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ABSTRACT 

Aims. We describe the observing strategy, data reduction tools and early results of a supernova (SN) search project, named SUDARE, 
conducted with the ESO VST telescope aimed at measuring the rate of the different types of SNe in the redshift range 0.2 < z < 0.8. 
Methods. The search was performed in two of the best-studied extragalactic fields, CDES and COSMOS, for which a wealth of 
ancillary data are available in the literature or public archives. We developed a pipeline for the data reduction and rapid identification 
of transients. As a result of the frequent monitoring of the two selected fields we obtained light curve and colour information for 
the transients sources that were used for the selection and classification of SNe by means of a especially developed tool. Eor the 
accurate characterisations of the surveyed stellar population we exploit public data and our own observations to measure the galaxy 
photometric redshifts and rest frame colours. 

Results. We obtained a final sample of 117 SNe, most of which are SN la (57%) and the remaining core collapse events of which 
44% type II, 22% type Iln and 34% type Ib/c. In order to link the transients, we built a catalog of ~ 1.3 x 10^ galaxies in the redshift 
range 0 < z < 1 with a limiting magnitude Kab = 23.5 mag. We measured the SN rate per unit volume for SN la and core collapse 
SNe in different bin of redshifts. The values are consistent with other measurements from the literature. 

Conclusions. The dispersion of the rate measurements for SNe-Ia is comparable with the scatter of the theoretical tracks for single 
(SD) and double degenerate (DD) binary systems models, therefore the data do not allow to disentangle among the two different 
progenitor scenarios. However, we may notice that among the three tested models, SD and two flavours of DD, either with a steep 
(DDC) or a wide (DDW) delay time distribution, the SD gives a better fit across the whole redshift range whereas the DDC better 
matches the steep rise up to redshift ~ 1.2. The DDW appears instead less favoured. The core collapse SN rate is fully consistent, 
unlike recent claims, with the prediction based on recent estimates of the star formation history, and standard progenitor mass range. 

Key words. Stars: supemovae: general - Galaxies: star formation - Galaxies: stellar content - surveys 


1. Introduction 


The evolution of the SN rate with redshift provides the observa¬ 
tional link between the cosmic star formation history (SFH), the 
initial mass function and the stellar evolutionary scenarios lead¬ 
ing to the explosions. Until recent times, the available measure¬ 
ments were limited to the local Universe and to sparse, some¬ 


times conflicting, high redshift measurements(e.g. Dahlen et al. 
2012[ Maoz et al.|2014[ and reference therein). The new gener¬ 


ation of panoramic detectors, now available in many observato- 


* Based on observations made with ESO telescopes at the Paranal 
Observatory under programme ID 088.D-4006, 088.D-4007, 089.D- 
0244, 089.D-0248, 090.D-0078, 090.D-0079, 088.D-4013, 089.D- 
0250, 090.D-0081 


ries, has substantially improved the survey capabilities and as a 
consequence the number of SN searches and rate measurements. 
Most of the efforts were devoted to type la SNe whose progenitor 
scenario is still strongly debated but the interest for core collapse 
SNe (CC SN) is also growing. 


The notion that measurement of the evolution of the SN la 
rate with redshifts, in combination with measurement of the 
SFH, can be used to constrain the SN la progenitor scenarios 
was first illustrated by Madau et al.|(|199^ (see also |Sadat et al.| 
1998[ Ruiz-Fapuente & Canal|1998|l. Early measurements were 


puzzling, showing a very rapid raise of the SN la rate up to red¬ 
shift ~ 1 and then a decline at higher redshift (Dahlen et al. 2004[ 
Barris & Tonry|[2006 Dahlen et al.||200^ . This implied a long 


delay time from star formation to explosion for the SN la progen- 
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itors, that appears conflicting with the indications derived from 


the rate measurement in the galaxies of the local Universe (Man- 
nucci et al.|2005|l. Subsequent measurements did not confirm the 


early results, but the issue is still debated (e.g. [Kuznetsova et al. 


2008 


2012 


Rodney & Tonryj 2010| Smith et al. [2012 Perrett et al. 


[Graur et al.|2014| and reference therein). 


For what concerns CC SNe, the early measurements show 
that, as expected, the rates tracks the cosmic star history (|Botti- 


cella et al.|2008 1, though the scaling factor seems to be a factor 


2 smaller than what was expected from the SFH. A possible ex¬ 
planation is that many dim CC SNe, are missed by SN searches 
( jHoriuchi et al. 201 l| l and/or part of the missing SNe may be hid¬ 
den in the dusty nuclear regions of star-burst galaxies ( [Dahlenj 
et al.|201^ and reference therein). Flowever, it is fair to say that 


the significance of the claimed discrepancy is still poor. 

In addition, we note that in most cases, the cosmic SN rates 
were derived from surveys designed to identify un-reddened 
SN la for the cosmological distance ladder. In these cases the 
specific observing strategy and/or candidate selection criteria 
may introduce biases in the event statistics that can be difficult to 
properly account. As a consequence, the SN rates derived from 
these surveys may be inaccurate. 

In order to use the rate evolution to constrain the SN pro¬ 
genitor scenarios, the knowledge of the properties of the parent 
stellar population is of fundamental importance, so that the vol¬ 
ume of Universe searched for SNe needs to be characterized in 
terms of the galaxy distribution as a function of redshift, mass 
and star formation history. 

Based on these considerations, we conceived a new SN 
search (Supernova Diversity and Rate Evolution, SUDARE, |Bot-| 
jticella et al.||2013] ), with primary goal to measure SN rates at 
medium redshift, that is 0.2 < z < 0.8. To combine the require¬ 
ments of good statistics (> 200 events) and the availability of 
ancillary data for the surveyed fields, we planned for a four year 
project to monitor of two very well known extragalactic fields, 
the Chandra Deep Field South (CDFS) and the Cosmic Evolu¬ 
tion Survey (COSMOS) fields. Thanks to the long term com¬ 
mitment of many different observing programs, extended multi¬ 
band photometry is available for these fields. These data allow an 
accurate characterization of the galaxy sample, which is crucial 
to infer general properties of the SN progenitors. 

The present paper describes the SUDARE survey strategy, 
the procedures for the identification of transients, and the SN 
candidate selection and classification after the first two years of 
observations. Then, we discuss the definition of the galaxy sam¬ 
ple and the procedure to derive photometric redshifts. Einally, 
we estimate the SN rates per unit volume at different redshifts 
and compare them to published estimates. A detailed study of 
SN rates as a function of different galaxy parameters will be pre¬ 
sented in a companion paper (Botticella et al., in preparation, 
hereafter PII). 

Along the paper, we will adopt Ha - 70km s ' Mpc *, Gm = 
0.3 and Da = 0.7. Magnitudes are in the AB system. 


2. The survey 

The SUDARE SN survey is performed using the VET Survey 
Telescope (VST, |^paccioli & Schipani|201 l| l equipped with the 
OmegaCAM camera ( |Kuijken|201 Ij l, that started regular opera¬ 
tion in October 2011 at ESO Paranal (Chile). The VST has a 
primary mirror of 2.6 m and a f/5.5 modified Ritchey-Chretien 
optical layout designed to deliver a large, uniform focal plane. 
The camera is equipped with a mosaic of 8 x 4 CCDs, each 
with 4kx2k pixels, covering 1 square degree with a pixel scale 
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of 0.214"pix *, allowing for an optimal sampling of the PSP 
even in good seeing conditions. The thinned CD44-82 devices 
from E2V have the advantage of an excellent quantum efficiency 
in the blue bands but with the drawback that the i and z-bands 
suffer for significant fringing contamination. 

Most of the observing time at this facility is committed to 
ESO public survey^ but a fraction of the time is dedicated to 
Guaranteed Time Observations (GTO) made available to the 
telescope and instrument teams in reward of their investments 
in the construction and installation of the instruments. 

SUDARE is a four-year program and this paper is devoted 
to the analysis of first two observing seasons for VOICE-CDPS 
and one season for COSMOS. We are now completing the mon¬ 
itoring of both fields for the subsequent two seasons. 

The time allocated to our project for monitoring VOICE- 
CDPS was from the VST and OMEGAcam GTO. The observing 
strategy is to span 4 deg^ in four pointings, with one pointing for 
each observing season (August to January). Here we present data 
for two of these pointings. 

To extend the photometric coverage, we implemented a syn¬ 
ergy with the VOICE (VST Optical Imaging of the CDPS and 
ESI Pields) project (Covone et al. in preparation). VOICE is a 
GTO program that has the aim to secure deep optical counterpart 
to existing multi-band photometry of selected fields. The multi¬ 
band catalog will be used to study the mass assembly and star 
formation history in galaxies by combining accurate photomet¬ 
ric redshifts, stellar masses and weak lensing maps. 

The monitoring of the COSMOS held relies instead on a pro¬ 
posal submitted for ESO VST open time (RI. Pignata). Por this 
held we maintained the same pointing coordinates from one sea¬ 
son to the next. This allows us the detection of transients with 
very long time evolution, but, due to the limited area probed, the 
data may be prone to cosmic variance. 

The survey strategy consists of monitoring the selected fields 
every three days in the r-band, excluding only ±5 days around 
full moon. The exposure time is 30 min with the aim of reaching 
a magnitude limit of 25* in average sky conditions. Each obser¬ 
vation is split into five 6 min exposures with a dithering pattern 
designed to All the gaps between the detector chips (that range 
from 25 to 85 arcsec). Because of the dithering, the effective area 
covered by combining the exposures of a given held is 1.15 deg^, 
although with a reduced signal to noise ratio at the edge. 

With a more relaxed cadence (3-4 times per month), we also 
planned for g and i-band exposures. With these observations we 
can measure the colour of the transients that is essential for their 
photometric classification and to obtain an estimate of the ex¬ 
tinction along the line of sight. Adapting to the rules for the 
observing blocks in service mode at ESO, the planned observ¬ 
ing sequence is g-r, r, r-i with a three-day interval between each 
block. To ensure a good image quality we required a maximum 
seeing of 1.2 " (PWHM) at the beginning of the exposure. This, 
along with the obvious requirement of clear sky, implied that the 
actual epochs of observations often deviate from the ideal plan 
mainly because of unsuitable sky conditions. 

Table lists for each held, the pointing coordinates along 
with the held size, the observing season, the number of available 
epochs in the different bands and the range and median value 
of the seeing measured for the r-band exposures. The full log of 
observations is given in Tab. j^where, for each epoch, we list the 
seeing (in arcsec) and mso, the magnitude corresponding with a 
transient detection efficiency of 50% (cf. Sect.[^. 

' http://WWW.eso.org/public/teles-instr/ 
surveytelescopes/vst/surveys/ 
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Table 1. Field coordinates and compact log of observations 


field 

R.A. (2000.0) Dec. 
hh:mm:ss dd:mm:ss 

field size 
deg2 

observing 

season 

epochs 

r g i 

seeing 

range [median] 

VOICE-CDFSl 

03:33:34.506 

-27:34:10.78 

1.15 

Aug2012-Jan2013 

29 

7 

11 

0.51-1.44 [0.89] 

VOICE-CDFS2 

03:29:02.654 

-27:34:00.70 

1.15 

Oct2011-Jan2012 

23 

6 

4 

0.51-1.46 [0.82] 

COSMOS 

10:00:28.600 

-1-02:12:21.00 

1.15 

Dec2011-Apr2012 

28 

7 

7 

0.50-1.20 [0.84] 


Table 2. Log of observations. For each field we report the epoch of observation (civil date and MJD), the seeing (FWHM in arcsec and mso the 
magnitude corresponding with a transient detection efficiency of 50%.) 


VOICE-CDFSl 

VOICE-CDFS2 

COSMOS 

Epoch 

MJD 

seeing 

mso 

Epoch 

MJD 

seeing 

niso 

Epoch 

MJD 

seeing 

mso 

2012-08-05 

r 

56144.38 

1.31 

23.0 

2011-10-20 

r 

55854.15 

1.17 

23.2 

2011-12-18 

r 

55913.30 

0.65 

23.2 

2012-08-13 

56152.37 

0.69 

23.9 

2011-10-25 

55859.34 

0.56 

23.8 

2011-12-22 

55917.27 

0.92 

23.5 

2012-09-02 

56172.23 

1.02 

21.7 

2011-10-28 

55862.16 

0.92 

23.5 

2011-12-27 

55922.24 

1.03 

23.3 

2012-09-05 

56175.20 

1.28 

21.8 

2011-10-30 

55864.16 

1.06 

23.4 

2011-12-31 

55926.28 

1.14 

23.3 

2012-09-08 

56178.31 

1.00 

23.0 

2011-11-02 

55867.10 

0.78 

23.5 

2012-01-02 

55928.32 

0.63 

23.7 

2012-09-14 

56184.25 

0.55 

23.9 

2011-11-04 

55869.15 

0.62 

23.4 

2012-01-06 

55932.30 

0.58 

23.7 

2012-09-17 

56187.29 

1.06 

23.2 

2011-11-15 

55880.06 

0.61 

23.6 

2012-01-18 

55944.20 

0.63 

23.7 

2012-09-20 

56190.22 

0.87 

23.1 

2011-11-18 

55883.29 

0.90 

23.4 

2012-01-20 

55946.25 

0.87 

23.4 

2012-09-22 

56192.22 

0.89 

23.1 

2011-11-21 

55886.23 

0.68 

23.5 

2012-01-22 

55948.25 

0.77 

23.6 

2012-09-24 

56194.24 

1.44 

22.6 

2011-11-23 

55888.28 

0.90 

23.4 

2012-01-24 

55950.28 

0.68 

23.8 

2012-10-07 

56207.34 

0.93 

23.1 

2011-11-26 

55891.28 

0.64 

23.8 

2012-01-27 

55953.19 

0.92 

23.5 

2012-10-08 

56208.24 

0.93 

22.9 

2011-11-28 

55893.30 

1.04 

23.2 

2012-01-29 

55955.20 

0.86 

23.3 

2012-10-11 

56211.32 

0.92 

23.1 

2011-12-01 

55896.22 

0.82 

23.6 

2012-02-02 

55959.31 

0.88 

23.3 

2012-10-14 

56214.71 

1.07 

23.7 

2011-12-03 

55898.24 

0.52 

23.8 

2012-02-16 

55973.25 

0.50 

23.6 

2012-10-17 

56217.21 

0.92 

23.7 

2011-12-14 

55909.21 

0.88 

23.3 

2012-02-19 

55976.11 

0.97 

23.3 

2012-10-21 

56221.19 

0.51 

23.8 

2011-12-17 

55912.25 

0.88 

23.2 

2012-02-21 

55978.15 

0.77 

23.6 

2012-10-25 

56225.11 

0.86 

22.9 

2012-01-14 

55940.15 

0.77 

23.3 

2012-02-23 

55980.18 

0.74 

23.6 

2012-11-04 

56235.14 

0.67 

23.6 

2012-01-18 

55944.06 

0.57 

23.9 

2012-02-26 

55983.12 

0.84 

23.7 

2012-11-06 

56237.26 

0.83 

23.5 

2012-01-20 

55946.134 

1.00 

23.4 

2012-02-29 

55986.040 

0.89 

23.2 

2012-11-08 

56239.27 

0.88 

23.6 

2012-01-23 

55949.11 

0.59 

23.5 

2012-03-03 

55989.196 

0.93 

23.2 

2012-11-10 

56241.30 

0.76 

23.6 

2012-01-25 

55951.12 

0.90 

23.7 

2012-03-06 

55992.115 

0.80 

22.6 

2012-11-20 

56251.14 

0.78 

24.0 

2012-01-29 

55955.06 

0.67 

23.6 

2012-03-13 

55999.033 

0.68 

23.5 

2012-12-03 

56264.21 

0.71 

23.4 

2012-02-02 

55959.08 

1.46 

22.4 

2012-03-15 

56001.045 

1.11 

23.1 

2012-12-07 

56268.05 

0.81 

23.8 


g 



2012-03-17 

56003.05 

0.92 

23.3 

2012-12-13 

56274.06 

0.55 

23.8 

2011-11-02 

55867.12 

0.91 

23.4 

2012-05-08 

56056.000 

0.71 

23.4 

2012-12-20 

56281.13 

0.96 

23.3 

2011-11-21 

55886.25 

0.89 

23.5 

2012-05-11 

56058.10 

0.85 

23.3 

2013-01-03 

56295.10 

0.68 

23.8 

2011-12-01 

55896.25 

1.05 

23.7 

2012-05-17 

56064.03 

0.75 

23.5 

2013-01-06 

56298.13 

0.91 

23.7 

2012-01-18 

55944.09 

0.56 

23.6 

2012-05-24 

56071.07 

1.20 

22.8 

2013-01-10 

56302.12 

0.89 

23.0 

2012-01-25 

55951.15 

1.56 

23.2 


g 




g 



2012-12-08 

56269.26 

1.02 

23.4 

2011-12-27 

55922.26 

1.13 

23.5 

2012-09-20 

56190.24 

1.17 

23.6 


i 



2012-01-22 

55948.27 

1.11 

23.6 

2012-10-11 

56211.35 

1.13 

23.4 

2011-11-02 

55867.15 

0.59 

23.1 

2012-02-02 

55959.33 

0.88 

23.8 

2012-10-21 

56221.22 

0.57 

23.7 

2011-11-21 

55886.28 

1.12 

22.7 

2012-02-16 

55973.28 

0.61 

23.9 

2012-11-06 

56237.28 

0.91 

23.6 

2011-12-01 

55896.27 

0.92 

22.8 

2012-02-26 

55983.14 

1.04 

23.8 

2012-12-07 

56268.07 

0.82 

24.1 

2012-01-18 

55944.11 

0.69 

22.9 

2012-03-17 

56003.07 

1.04 

23.6 

2012-12-09 

56270.28 

1.08 

23.5 





2012-05-09 

56056.02 

0.78 

23.8 

2013-01-06 

56298.16 

1.14 

23.7 






i 




i 







2011-12-27 

55922.28 

0.93 

24.1 

2012-08-13 

56152.39 

0.54 

23.0 





2012-01-22 

55948.30 

0.95 

23.2 

2012-09-02 

56172.26 

0.98 

21.1 





2012-02-02 

55959.36 

1.05 

22.7 

2012-09-08 

56178.34 

1.45 

22.0 





2012-02-16 

55973.30 

0.60 

23.3 

2012-09-17 

56187.31 

0.91 

22.7 





2012-02-26 

55983.17 

0.88 

24.1 

2012-09-24 

56194.27 

1.42 

22.2 





2012-03-17 

56003.10 

0.53 

24.2 

2012-10-15 

56215.19 

0.89 

22.9 





2012-05-09 

56056.05 

0.62 

23.4 

2012-10-25 

56225.14 

0.85 

22.2 









2012-11-08 

56239.29 

0.84 

22.8 









2012-11-20 

56251.16 

0.69 

23.2 









2012-12-20 

56281.15 

0.79 

22.7 









2013-01-10 

56302.15 

0.97 

22.0 
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The SN search was complemented by 3 runs of one night 
each at the ESO-VLT for the spectroscopic classification of a 
dozen candidates. These observations, described in Sect. |4.3[ 
were intended as spot checks of the SN photometric classifica¬ 
tion tool (Sect. 0- 

We remark that, as by-product, the SUDARE data archive 
was also used to explore the performances and completeness of 
AGN detection via variability (|De Cicco et al.||2015[ |Falocco| 
[etal]2^. 


2.1. Image calibration 


The raw data were retrieved from the ESO archive and trans¬ 
ferred to the VST data reduction node in OAC-Naples. Here the 
first part of the data reduction is performed using the VSTTube 
pipeline ( |Grado et al.|2012| l. A description of the VSTtube data 
reduction process is reported in De Cicco et al. ( 2015| l. In short, 
the pipeline first performs flat fielding, gain harmonisation and 
illumination correction and all images for a given field are reg¬ 
istered to the same spatial grid and photometric scale. Finally 
the dithered images for one epoch are median averaged produc¬ 
ing one stacked image. The pipeline also delivers weight pixel 
masks tracking for each pixel the number of dithered exposures 
contributing to the combined image after accounting for CCD 
gaps, bad pixels and cosmic rays rejection. 

The pipeline was also used to produce deep stacked images 
by combining all the exposures in a given filter with the best 
image quality, those with seeing < 0.8". These stacked images, 
reaching a limiting magnitude ~ Imag fainter than good sin¬ 
gle epoch exposures (the 3-cr magnitude limits are 26.2, 25.6, 
24.9 mag for r,g and i-band respectively) were used to extract 
galaxy photometry to complement the public multi-band cata¬ 
logs. (Sect.|^. 


2.2. Transient detection 


For the detection of transient sources and the selection of SN 
candidates, the mosaic images were processed with an ad-hoc 
pipeline. This is a collection of python scripts that makes use 
of pyraf and tasks and incorporate other publicly avail¬ 

able software for specific tasks, in particular SExtractor ( |Bertin| 
|& Arnouts|1996|l for source detection and characterisati on, hot- 


pant. 


for PSF match and image difference, daophot (Stetson 


|1987 ]( for accurate point spread function (PSF) fit photometry, 
iti/f^for catalog handling and mysql for the transient database. 
The flowchart of the SUDARE pipeline is the following: 


1. We produced a mask for saturated stars that is combined with 
the weight map produced by VSTTube, to build a bad pixel 
mask for each mosaic image. Those flagged as bad pixels 
were excluded from further analysis. 

2. For each image we computed the difference from a selected 
template image. This required the derivation of the convolu¬ 
tion kernel matching the PSF of the two images. The method 
is described in |Alard| ( |2000) l though we used the hotpants 
implementation. We note that for the first observing season 
we do not have earlier templates and we therefore used as 
templates, images acquired on purpose a few months after 
completion of the transient survey campaign. 

^ http: //WWW. stsci.edu/institute/software_hardware/ 

^ A package provided by A. Becker (http://www.astro. 

Washington.edu/users/becker/v2.8/hotpants.htmlI 

** http: //WWW. star. hr is. ac. uk/~mht/stilts/ 
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3. The transient candidates were identified, using SExtractor, 
as positive sources in the difference image. Depending on 
the image quality and detection threshold, the candidate list 
starts with several thousands objects for each epoch. Most 
detections are artefacts due to poorly masked CCD defects, 
poorly removed cosmic rays, residual from the subtraction of 
bright sources, reflection ghosts from bright sources, etc. 

4. The transient candidates were ranked based on a custom 
algorithm that uses a number of measured sextractor met¬ 
rics of the detected sources. The most informative param¬ 
eters and the ranking scores were selected and calibrated 
through extensive artificial star experiments. In these exper¬ 
iments, a number of fake stars were placed in the search im¬ 
age that is then processed through the detection and ranking 
pipeline. The success rate of artificial star recovery was com¬ 
pared to the number of residual spurious sources. We found 
that the most informative parameters are the source FWHM, 
flux_ratio, isoarea and magnitudes, measured at different 
apertures. With a proper selection of these parameters, we 
can drastically reduce the number of spurious events while 
limiting the number of good candidates improperly rejected. 
The performance of the ranking algorithm depends on the 
image quality. On average, we found that we can eliminate 
~ 95% of the spurious transients at the cost of losing ~ 5% 
of good candidates. Correction for the lost SNe is incorpo¬ 
rated in the detection efficiency (c.f. Sec. since we used 
the same algorithm to select real and fake SNe. After this 
selection, typically a hundred candidates per field and epoch 
are left to the next step of human inspection and validation. 

5. To associate each transient to its possible host galaxy, we 
cross correlated the transient list with the galaxy catalogs 
derived from the deep r-band stacked image (cf. Sec. [^. A 
galaxy is adopted as the host for a given transient when the 
latter appears engulfed in the galaxy boundaries. The bound¬ 
aries are those of the ellipse defined by the SExtractor pa¬ 
rameters CXX, CYY and CXY through the equation: 

CXX(.r - Xcf + CYY(y - y,f + CXY{x - x^Xy - ye) = 

where Xc, yc is the galaxy center and following the SExtrac- 
toEs manual, we assume that the isophotal limit corresponds 
to R = 3. 

In four cases the transient/host galaxy pairing was ambigu¬ 
ous because of the overlap of the ellipses of different galax¬ 
ies. In these cases we also considered the consistency of the 
host galaxy redshift with the indication of SN photometric 
classification (c.f. Sect.|^. 

Only for two transients no counterparts was detected in the 
deep stacked image (cf. Tab.|^. 

6. The available information for the best ranked candidates are 
posted on WEB pages where the user can inspect the images 
and candidate metrics for the search and template epochs 
and selects the good candidates assigning each of them a 
preliminary classification according to different classes (SN, 
AGN, variable star, moving object). The selected candidates 
are then archived in a mysql database. 

7. For all selected candidates, we derive accurate light curves 
measuring the source magnitude at all available epochs. We 
measure both aperture and PSF-fit photometry on the origi¬ 
nal search images and on the difference images. We verified 
that the PSF-fit gives more reliable measurements that plain 
aperture integration mainly because PSF photometry is less 
sensitive to the background noise. 
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The transient search process was performed in the r-band for 
each epoch. Because of the dense temporal sampling, in gen¬ 
eral good SN candidates will have multiple detections in the 
database. In principle, we can easily implement a candidate se¬ 
lection based on multiple occurrences of a given source in the 
database that would further reduce spurious candidates. How¬ 
ever, at the present stage of the project we adopted a conserva¬ 
tive approach accepting the burden of the visual examination of 
many candidates to maximise the completeness fraction. 

3. Detection efficiency 

To derive the SN rates we need to obtain an accurate estimate 
of the completeness of our search. This is done by extensive ar¬ 
tificial stars experiments exploring a range of magnitudes and 
positions on the images. 

For every search image, we first obtain the PSF from the 
analysis of isolated field stars. Then, a number of fake stars of 
a given magnitude, generated scaling the PSF and adding the 
proper Poisson noise, were injected on the search image. To 
mimic the range of properties of real sources, three different cri¬ 
teria were adopted for positioning the fake stars, with roughly 
the same number of star for each class. The three classes are: 



Fig. 1. Transient detection efficiency as a function of magnitude for 
the r-band observation of COSMOS on 2012/03/15. The dots are the 
averages of three artificial star experiments with the errorbar being the 
dispersion whereas the line is the adopted efficiency curve after the fit 
with the = 93%, mag^o = 23.0 mag and fi = 6.1). 


- events associated with galaxies. From the source catalog on 
the field (cf. Sect.|^ we picked up a random sample of galax¬ 
ies and one fake star was placed in each of them. The position 
inside the galaxy was chosen randomly following the distri¬ 
bution of r-band flux intensity. 

- events coincident with persistent, point-like sources. Fake 
stars were added at the same position of existing sources in 
the field. This mimics SNe in the nucleus of compact host 
galaxies, variable AGNs and variable stars. 

- events with no counterpart in the template image. These were 
placed at random positions across the field of view irrespec¬ 
tive of existing sources or to the possible coincidence with 
CCD defects or gaps. 

The images with fake stars were processed with the search 
pipeline and the number of detected events surviving after the 
ranking procedure were counted. The fraction of detected over 
injected events gives the detection efficiency for the given mag¬ 
nitude. The experiment is repeated sampling the magnitude 
range of interest (18 < r < 26) to derive the detection efficiency 
as a function of magnitude. 

We inject 500 fake stars per experiment per image and repeat 
the experiment 5 times for a given magnitude. In fact, adding 
a large number of fake stars could bias the computation of the 
convolution kernel, hence the image subtraction process, making 
the experiments less reliable. 

An example of the derived detection efficiencies as a function 
of magnitude for one epoch and field is shown in Fig. where 
the errorbars show the dispersions out of three experiments. As 
can be seen in Fig. the detection efficiency as a function of 
magnitude can be represented by the following analytical func¬ 
tion: 


rvc r .17 /'arctan^(ma 05 o-mag) 1 \ 

DE = DE,„ax I - - -+ - I (1) 

where DE^ax is the maximum value of the detection effi¬ 
ciency, magso is the magnitude corresponding with the 50% drop 
of the detection efficiency and fi measures the decline rate of the 


DE. The best fit parameters are determined through least square 
minimisation of the residuals (e.g. Fig. [^1. Note that the maxi¬ 
mum detection efficiency is ~ 95%, even at the bright magnitude 
end. This is because not all the pixels of the image are useful and 
indeed the bad pixel mask (cf. step 1 of the SUDARE pipeline) 
flags 5 - 10% of the image area. 

We verified that the detection efficiencies measured indepen¬ 
dently for each of the three classes of fake stars described above 
are similar with a dispersion of mag^Q values of ~ 0.3 mag and 
hereafter we will use the average values. Also, we found that the 
values of DE^ax and /? are very similar for each epoch and field 
with a mean value of 95% + 3 and 4 + 2, respectively. 

The artificial star experiment described above was repeated 
for all epochs and fields of the search, and the resulting detection 
efficiencies are used in the rate calculation. 


4. Transient classification 


The result of the transient search was a list of ~ 350 SN candi¬ 
dates. A fraction of these transients are coincident with persis¬ 
tent sources and therefore can be SNe in the nucleus of the host 
galaxies but also variable stars or AGNs. In fact, after the anal¬ 
ysis of |Falocco et al. ( 2015| l and De Cicco et al. (20151, three 
candidates with slow evolving light curves were found to co¬ 
incides with X-ray sources. These were classified as AGN and 
removed from the SN candidate sample. We also removed from 
the SN candidate list all transients with spectroscopic or photo¬ 
metric redshift z > 1 (~ 25%). 

We used the measured light curve and colour evolution to 
constrain the nature of the transients and to classify SNe in their 
different types. Conservatively, only candidates with at least 5 
photometric measurements at different epochs (even in different 
filters) were considered. Because of the frequent monitoring of 
our survey, this criterium excludes only four candidates with a 
negligible effect on the SN counts. 

The photometric classification is more reliable if the redshift 
of the host galaxy is available. When the transient is not associ¬ 
ated with a host galaxy or when the host redshift is not available. 
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the redshift is left as a free parameter in the transient light curve 
fitting (see next section). 

4.1. Photometric classification of SNe 

For the photometric classification we used a tool developed 
for the SUDARE project. The tool compares the SN candidate 
multi-color light curves with those of SN templates, and iden¬ 
tifies the best-matching template, redshift, extinction and lumi¬ 
nosity class. The tools was developed following the strategy of 
the SN classification tool PSNID ( |Sako et ar]|201 l| l. We devel¬ 
oped our own tool because we want to explore different priors 
for the fitting parameters and different classification scheme. 

We collected a sample of templates for different SN types 
for which both multicolour light curves and sequence of spectra 
are available (Tab.[^. The spectra are needed to estimate the K- 
correction. The templates were retrieved from a database of SN 
light curves and spectra that we have collected in the study of 
SNe at ESO and ASIAGO Observatow (the template spectra can 
be downloaded also from WISEREl^ the SN spectra database; 
|Yaron & Gal-Yam|2012| l. 

The templates were selected to represent the well established 
SN types, namely la, Ib, Ic, Ilb, II Plateau and Linear, Iln, with 
the addition of representative peculiar events (see individual ref¬ 
erences for details). In particular, we included SN 2008es as rep¬ 
resentative of the recently discovered class of very luminous SNe 


(SLSN, Quimby et al. 2007[ 

20131 Gal-Yam 

2012 1 that, although 

intrinsically very rare, may 

)e detected in high redshift searches 


because of their large volume sampling. 

The steps of the photometric typing are; 

- for each template, we derived /T-correction tables as a func¬ 
tion of phase from maximum and redshift (in the range 
0 < z < 1). /T-corrections are obtained as the difference of 
the synthetic photometry measured on the rest frame spec¬ 
tra and on the same spectra properly redshifted. The redshift 
range for which we can derive accurate /T-correction is lim¬ 
ited by the lack of UV coverage from most templates. This 
is a problem in particular for the g-band where we are forced 
to accept uncertain extrapolations. 

- the K-corrected light curves of template SNe were used to 
predict the observer frame light curves in the gri-bands ex¬ 
ploring the 0 < z < 1 redshifts range and the -0.3 < Eb-v < 
1 mag extinction range (the negative lower limit for the Eg-v 
range allows for uncertainties in the template extinction cor¬ 
rection and for variance in the intrinsic SN colour). With the 
goal to minimise the uncertainties on the /T-correction, the 
template input band was taken to best match the observer 
frame band for the given redshift, e.g. we use the template 
V, B, U bands to predict the observer frame r-band light curve 
of SNe at redshift z ~ 0.1, 0.4, 0.7, respectively. 

- we estimate the goodness of the fit of the template to the 
observed light curve computing the sum of the square of flux 
residuals weighted by the photometric errors (x^) for each 
simulated light curve of the grid. Besides the redshift and 
extinction ranges we explore a range of epochs of maximum, 
Tmax (the initial guess is the epoch of the observed r-band 
brightest point) and of intrinsic luminosities, A(jj.) (allowing 
for a +0.3 mag flux scaling of the template). The residuals 
for all bands are summed together and therefore each band 
contributes to the overall;^^^ with a weight proportional to the 
number of measurements. 

^ http://wiserep.weizmann.ac.il 
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for the selection of the best fitting template we use Bayesian 
model selection (e.g. Poznanski et al.|2007a| 


Kuznetsova & 


Connolly 2007 [ Rodney & Tonry 2009 | l. In particular fol¬ 
lowing Sako et al. ( |2011| l we compute for each SN type the 
Bayesian evidence: 


E 


type - 


^ I P(z)e ^^^^dzdAvdT^^dAfjj.) 


where the fitting parameters are the redshift z, with P{z) its 
probability distribution, the extinction Ay, the time of maxi¬ 
mum Tmax and the flux scaling factor A{p). 

The spectroscopic redshift is used as a prior if available and 
in this case for P{z) we adopt a normal distribution centered 
at the spectroscopic redshift and with cr - 0.005. Other¬ 
wise, if a photometric redshift estimate is available, we use 
as redshift prior the P(z) provided by the photometric red¬ 
shift code (cf. Sect. |6.3.1| l. In the worst case, either when 
the host galaxy is not detected, or when the photometric red¬ 
shift is poorly constrained, we adopt a flat prior in the range 
0 < z < 1. 

In all cases we adopted flat prior for the extinction distribu¬ 
tion and for the flux scaling. 


More critical is choice of templates. As emphasised by Rod- 
|ney & Tonry| ( |200^ , the Bayesian approach relies on an ap¬ 
propriate template list that should be as complete as possible 
but, at the same time, avoid duplicates. When the template 
list includes rare, peculiar events, especially if they mimic 
the properties of a more frequent SN type, it is appropriate 
to use frequency priors. Alternatively, for specific applica¬ 
tions one may exclude ambiguous cases or rare, peculiar SN 


types (c.f.|Sako et al.|201 l|l from the template list. 

Our template list, given in Tab. is intended to represent 
the full range of the most frequent SN type with a number of 
templates for each class that is broadly consistent with their 
frequency in a volume limited SN sample ( |Li et al.|201 lb| l. 
After that we adopted flat priors for the relative rate of each 
template within a given class and for the relative rates of the 
different SN types. 

We computed the Bayesian probability for each of the main 
SN types as: 


Ftype — 


-type 


Ela + Eib/c + Eli + (Eslsn) 


Notice that for the purpose to assign probability, we merged 
regular type II and type Iln templates. However, in the sub¬ 
sequent analysis we indicated when the best fitting template 
(the one with the highest probability) is a type Iln. Also, after 
verifying that none of our candidate has a significant proba¬ 
bility of matching a SLSN, the corresponding template was 


dropped from the fitting list and the Eslsn term in Eq. 4.1 


cancelled. This was done to allow a direct comparison with 
SNANA (see next Section). 

for the most probable SN type, we record the best fitting tem¬ 
plate along with the fit parameters corresponding to the 
minimum. Of the about 250 transients, 117 were classified 
as SNe. Most of the remaining have erratic light curves con¬ 
sistent with that of AGNs. 

An example of the output of the SN typing procedure is 
shown in Eig. 


Table lists our SNe; for each event we report the coordi¬ 
nates, the host galaxy redshift if available (col. 4), for photomet¬ 
ric redshift, the 95% lower and upper limits of the P(z) (col. 5), 
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Fig. 2. Example of the output of the SN typing procedure. The top panel 
shows the observed r-band light curve that in this case is compared with 
the template A'-corrected 6-band light curve. The bottom panels show 
the observed light curve and template fit for the the g'-band (left) and 
;-band (right). Blue dots are the SN candidate observed magnitudes (ar¬ 
rows indicate upper limits) while red open circles are the template pho¬ 
tometry. The legend identifies the best fitting template and parameters. 


cosmos r 20111218 20120508 lA 



0 50 100 0 50 100 
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the method of redshift measurement (col. 6) the most probable 
SN type (col. 7) and corresponding Bayesian probability (col 8), 
the best fitting template (col. 9), redshift (col. 10), extinction 
(col. 11), flux scaling (col 12) and epoch of maximum (col. 13). 
We also list the (col. 14), the number of photometric mea¬ 
surements with, in parenthesis, the number of measurement with 
SNR > 2 (col. 15) and the integrated right tail probability of the 
distribution (B^j, col 16). 


In a number of cases the P^i probability is fairly low (15 SNe 
have Px^ < 10“"^). In some cases this is due to one/two deviant 
measurements, in other cases there is evidence of some variance 
in the light curve not fully represented by the adopted template 
selection. We have to consider the possibility that these events 
are not SNe. 


Also for some candidates with a small x^, the number of ac¬ 
tual detections (photometric measurements with signal to noise 
ratio S NR > 2) is also small that it is not possible to definitely 
assess the SN nature of the transient (for 12 candidates the num¬ 
ber of detection is Ndet <-2)- 


To these probable SNe (indicated with PSN in the last col¬ 
umn of Tab. we will attribute a weight 0.5 in the rate calcula¬ 
tion. The impact of the arbitrary thresholds for Px^ and the Adet 


and the adopted PSN weight will be estimated in Sect. 7.3 


To evaluate the uncertainties of our classification tool in the 
next two sections we compare our derived SN types with a) 
photometric classifications obtained using the public software 
package SNANA ( Kessler et al.||200'9) l and b) with the spectro¬ 
scopic observations of a small sample of "live" transients, ob¬ 
served while still in a bright state. 


Table 5. Comparison of photometric classification with the different 
tools. In parenthesis we report the events labelled as probable SNe. Eor 
column 2 and 3 we count the SN with respect to the most probable SN 
type. Column 4 lists instead sum the Bayesian probability for each SN 
type (in parenthesis is the number of PSN). 


SUDARE 

SNANA 

Bayesian 

la 

67 

72 

64.7 

(12.7) 

II 

22 

21 

23.2 

(6.5) 

Iln 

11 

7 

10.7 

(2.9) 

Ib/c 

17 

15 

18.6 

(4.8) 

All 

117 

115 

117.0 

(27.0) 


4.2. Comparison with the photometric dassifications by 
PSNID in SNANA 

With the aim to check our procedure and evaluate the related 
uncertainties, we performed the photometric classification of our 
SN candidates using the public code PSNID in the 5A4A4|^ 

Also for the fit with PSNID we set the host galaxy redshift as 
a prior with the same range of uncertainty as in our procedure. 
In this case however, also for photometric redshift we assume a 
normal distribution for P(z) with the cr provided by the photo¬ 
metric redshift code. 

A comparison of the classifications obtained with the two 
tools is illustrated in Fig. The pie chart shows the SN clas¬ 
sifications in the four main types using the SUDARE tool and 
the sectors with different colours within a given wedge show the 
PSNID classifications. For a two events, marked in grey, the fit 
with PSNID fails. 

The figure shows that the identification of SN la is quite con¬ 
sistent (92% of the type la classified by our tool are confirmed 
by PSNID) and there is a good agreement also for the normal 
type II (77% of the classifications are matched). The agreement 
is poor for type Ib/c and for type Iln (only 40% are matched 
in both cases). The latter result is not surprising considering the 
wide range of luminosity evolution: for these classes of SNe the 
choice of input templates is crucial. 

However, we remark that despite the discrepancy in the clas¬ 
sification of individual events of specific sub-types, there is an 
excellent agreement of the event counts in each class, but for 
type Iln, as shown in Tab. This implies that, as far as the SN 
rates are concerned, using either classification tools makes a lit¬ 
tle difference, with the exception the exception of type Iln where 
the difference is ~ 40%. 

In Tab|^we also report the SN count for the different classes 
using the Bayesian probability. It appears that with respect to the 
count of the most probable type, the number of SN la is slightly 
reduced while the number of Ib/c increases. This is not unex¬ 
pected given the similarity of the light curves of type la and Ib/c 
(in many cases an event can have a significant probability of be¬ 
ing either a type la or a type Ib/c) and the fact that type la are in¬ 
trinsically more frequent than Ib/c. The effect however is small, 
< 5% in both cases. 



size@update\enc@updateSNANA}versionlO_39k. mplementation 
( |Sako et al.|201 l[[Kessler et al.|2009| l. Overall, the approach of PSNID 
is similar to that adopted here; besides the implementation of the 
computation algorithm, the main difference is in the template list. In 
particular, for SN la we adopted the fitting set-up of |Sako et al.| ( |20lT) l, 
while for CC SNe we used the extended list of 24 templates available 
in the SNANA distributioij^ 
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Table 3. List of template SNe used for the SUDARE’s SN photometric classification tool 


SN type vel m- M reference 

[kms-'] 


1990N 

la 

998 

31.73 

Lira et al. 

(1998i; 

Mazzali et al. 

1993D 



1992A 

la 

1845 

31.14 

Krrshner et al. ( !: 

93 i; Altavilla et al. 

(^2004 

n 


1994D 

la 

450 

30.92 

Richmond et al. d 

1995 T; Patat et al. ( 

996: 



2002bo 

la 

1289 

31.77 

Benetti et al. (|20ti 

ITIi 





1999ee 

la bright 

3407 

33.42 

Stritzmger et al. d 

ZDO 2 I 1 Hamuy et al.l 

2 OO 2 I 1 


1991T 

la bright 

1732 

30.74 

Lira et al. 

rrwsi; 

Altavilla et al. 

(poo 

4li; Ruiz-Lapuente et al. 

I 992 I 1 

1991bg 

la faint 

913 

31.44 

Filippenko et al. ( 

1992 1 ; Leibundgut 

;t aLl 

T993|; 

liiratto et al 

IT996I 

2000cx 

la pec 

2421 

32.39 

Candia et al. (200 

3 1 ; Matheson et al. 

2008 

1 


2002cx 

la pec 

7183 

35.09 

lA et al. 

(1200311 






1987A 

11 

320 

18.48 

Gatchpole et al. ( 

989 1 ; Hamuy et al. 

(1990 1 ; ESO/Asiago Archive 

1992H 

IIL 

1021 

30.97 

Clocchiatti et al. i 

1996a 1 





2009bw 

IIP 

1155 

31.45 

Inserra et al. dZOlZll 





1999em 

IIP 

710 

29.47 

Hamuy et al. (20t 

rr 1 ; Leonard et al. (2 

002a 

1 ; Elmhamdi et al. ( 

20031 

2004et 

IIP 

40 

28.85 

Maguire et al. (2(j 

TOi 





1999br 

II faint 

1021 

30.97 

Pastorello et al. 

10041 





1999gi 

II faint 

592 

29.80 

Leonard et al. (|2U 






2005cs 

II faint 

600 

29.26 

Pastorello et al. (jl- 

!006 2009|l 





1993J 

lib 

-35 

27.80 

Filippenko et al. ( 

1992 1 ; Richmond et 

al. (1994i; 

Barbon et al. 

(19951 

2008ax 

lib 

579 

29.92 

Pastorello et al. (1^ 

!008 1 ; laubenberger et al. 

(|2U11 

1 

1997cy 

n pec 

17700 

37.03 

Turatto et al. (20L 

iOi 





1998S 

Iln 

895 

31.18 

Fassia etal. (2001 






2010jl 

Iln 

3207 

34.92 

Pozzo et al. (2004 






2005gj 

Iln (la) 

17988 

37.15 

Aldermg et al. (^20O6 \ 





2008es 

SFSN-II 

0.202* 

39.70 

Gezari et 

sr^oiT 






2009jf 

Ib 

2379 

32.65 

Valenti et al. (201 

Ti 





2008D 

Ib (XRF) 

1955 

32.29 

Mazzali et al. (20 

OSi 





19941 

Ic 

461 

29.60 

Wheeler et al. (|iy 

94 1 ; Clocchiatti et a 

lj( 1996b 1 


1998bw 

Ic (GRB) 

2550 

32.76 

Galama et al. (19' 

:^ 81 ; Patat et al. 

(poo 

1 



2004aw 

Ic 

4742 

34.17 

laubenberger et a 

1(20051 - 





2007gr 

Ic 

492 

29.84 

Valenti et al. (|200 

8 Hunter et al. (20 

09 1 




* redshift instead of velocity 



iF-jj unkw 

Fig. 3. Comparison of the SN classifications obtained with the different 
tools. The exploded wedges are the SN type fractions obtained with our 
SUDARE tool and the coloured sectors are the SNANA classification. 


4.3. Comparison with spectroscopic ciassification 

For a small sample of the SN candidates we obtained immediate 
spectroscopic classification. 

Observations were scheduled at the ESO VLT telescope 
equipped with FORS2 at three epochs for a total allocation of 
2 nights. The telescope time allocation, that was fixed several 
months in advance of the actual observations, dictated the choice 
of the candidates; we selected transients that were "live" (above 
the detection threshold) at the time of observations and among 
these we gave a higher priority to the brightest candidates with 
the aim to secure a higher S/N for the spectra. 

For the instrument set-up we used two different grisms, 
namely GRJS_300V and GRJS_300I covering the wavelength 
range 400-900 nm and 600-1000 nm respectively with similar 
resolution of about 1 nm. The choice of the grism for a particu¬ 
lar target was based on the estimated redshift of the host galaxy 
with the GR1SM_300I used for redshift z > 0.4. 

We were able to take the spectrum of 17 candidates. Spectra 
were reduced using standard recipes in IRAF. In three cases the 
S/N was too low for a conclusive transient classification and we 
were only able to obtain the host galaxy redshifts. Four of the 
candidates turned out to be variable AGNs, in particular Seyfert 
galaxies at redshifts between 0.25 < z < 0.5. We stress that, 
to maximise the chance of obtaining useful spectra, we tried to 
observe the candidate shortly after discovery. This means that at 
the time of observations we had not yet a full light curve and 
hence a reliable photometric classification. Eventually, all the 
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Fig. 4. The observed spectrum of SN 2012gs obtained with FORS2 on 
MJD 56252.0 (black line) is compared with that of the SN Ic 2007gr at 
phase -9d (top) and of the type la SN 199IT at phase +14d (bottom). 
In both cases it is adopted for SN 2012gs a redshift z = 0.5 as measured 
from the narrow emission lines of the host galaxy. 


four AGN exhibit a erratic luminosity evolution that, if known 
at the time of spectroscopic observations, would have allowed 
us to reject them as SN candidates. 

Ten transients were confirmed as SNe and their spectral type 
were assigned through cross-correlation with libraries of SN 
template spectra using GEneric cLAssification TOol (GELATO, 
[Harutyunyan et al.|2008|l and t he Supernova Identification code 
(SNID, |Blondin & Tonry|2007[ l. The spectroscopically classified 
SNe, identified with a label in Tab. turned out as 6 type la, 2 
type Ic, 1 type II and a type Iln. In all cases the SN type was 
coincident with the independent photometric classification with 
one exception (SN 2012gs) classified Ic from spectroscopy and 
la from photometry. 

As shown in Fig. the spectrum of SN 2012gs can be fit¬ 
ted both by a template of type Ic SN well before maximum or 
by a type la SN two weeks after maximum, in both cases the 
redshift was z ~ 0.5. On the other hand, when we consider the 
light curve (Fig.[^ it results that the spectrum was obtained two 
weeks after maximum and therefore the first alternative can now 
be rejected. Therefore, revising the original spectroscopic clas¬ 
sification, SN 2012gs is classified as type la. 


4.4. Classification uncertainties 


The comparison of the photometric and spectroscopic classifica¬ 
tions, even if for a very small sample, confirms that photometric 
typing is a reliable in particular when the redshift of the host 
galaxy is known. For our photometric tools we did not yet per¬ 
formed a detailed testing that instead has been performed for 
PSNID. In particular, |Sako et ^ ( |2011 1 show that PSIND can 
identify SN la with a purity of 90%. This appears consistent 
with the results obtained from the comparison of PSNID and 
SUDARE tools. More difficult is to quantify the performances 
of photometric classification for CC SN both because of the lack 
of suitable spectroscopic sample s (|Sako et al.|2011|l and for the 
limitation of simulated samples ( [Kessler et al.||201^ . From the 
comparison of the CC SN classification of the PSNID and SU¬ 
DARE tools we found differences in the individual classifications 
of 25% for type II events and 40% for type Ib/c events. These 


Fig. 5. SN 2012gs light curve fit obtained using our tool. The best match 
is obtained with SN 199IT and the maximum is estimated to occur on 
MJD 56235.1. 



discovery mag 


Fig. 6. Histogram of the r-band SN magnitudes at discovery. The dotted 
line shows the distribution of mjo (the magnitude where the detection 
efficiency is 0.5) for the r-band observations (as reported in Tab 2). 


should be considered as lower limit of the uncertainty because 
the two codes adopt similar approaches the main difference be¬ 
ing the choice of templates. On the other hand, the discrepancy 
on the overall SN counts of a given type is much lower, typically 
a few percents, though for type Iln it is about 40%. Based on 
these considerations and waiting for a more detailed testing we 
adopt the following uncertainties for SN classification: 10% for 
la, 25% for II, 40% for Ib/c and Iln. 


5. The SN sample 

As a result of the selection and classification process, we ob¬ 
tained a sample of 117 SNe, 27 of which are marked as probable 
SNe (PSNe). The distribution of the SN apparent magnitude at 
discovery, plotted in Fig.[^ shows a peak at mag r = 23.5 - 24 
that is consistent with what expected given the detection efficien¬ 
cies (see Sect|^. 

We found that 57% of the SNe are of type la, 19% of type 
II, 9% of type Iln and of 15% type Ib/c. We notice that the per- 
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centage in different subtypes is quite close to the fraction of SN 
types in magnitude limited samples. For instance the updated 
Asiago SN Catalog] includes 56% la, 27% II, 4% Iln and 10% 
Ib/c (counting only SNe discovered since 2000) with some dif¬ 
ference from our sample only for the most uncertain events clas¬ 
sified as type Iln. The result is encouraging when considering 
that we did not make any assumption on the fraction of the dif¬ 
ferent SN types in our typing procedure. This also implies that 
the relative rates of the different SN types are similar in the local 
Universe and at z ~ 0.5. 

At the same time the SN type distribution in our sample is 
very different from that derived in a volume limited sample like 
that derived for the LOSS survey ( |Li et al. |201 lb[ ) which gave 
the following SN type fraction; la 24%, Ib/c 19%, II 52 % and 
Iln 5%. The much higher fraction of SN la in our sample is a 
natural consequence of the high luminosity of SN la compared 
with other types that makes it possible to discover SN la in a 
much larger volume. This also explains the SN redshift distribu¬ 
tion shown in Fig. |7] While SN la are found to z ~ 0.8, the red¬ 
shift limit for the discovered SN II is only z ~ 0.4. The relatively 
rare but bright type Iln are, on average, discovered at higher red- 
shifts. 


6.1. COSMOS field 


A photometric catalog of the sources in the COSMOS field has 
been produced by |Muzzin et al.J _( |20T3] l, and available trough 
the Ultra VISTA survey Web sit^ The catalog covers an area 
of 1.62 deg^ and encompasses the entire 1.15 deg^ area moni¬ 
tored by SUDARE. The catalog includes photometry in 30 bands 
obtained from: /) optical imaging from Subaru/SuprimeCam 
(grizBV plus 12 medium/narrow bands IA427 - IA827) and 
CFHT/MegaCam (u*) (|Taniguchi et al.|2007l|Capak et al.|20071l, 
ii) NIR data from VISTA/ VIRCAM (YJHK hands, |McCracken| 
et al. 2012|l, Hi) UV imaging from GALEX (FUV and NUV 


channels, IMartin et al.|2005all, iv) MIR/FIR data from Spitzer’s 
IRACh-MIPS cameras (3.6, 4.5, 5 .8, 8.0, 24 and 70, 1 60 pm 
channels from Sanders et al.|2007 and Frayer et al.|2009) l. 

The optical and NIR imaging for COSMOS have compara¬ 
ble though not identical PSF widths (FWHMs are in the range 
0.5" - 1.2"). For an accurate measurement of galaxy colours, 
|Muzzin et al.| ( |2013] l performed the PSF homogenisation by de¬ 
grading the image quality of all bands to the image quality of 
the band with worst seeing (with seeing of 1"- 1.2"). Source 
detection and photometric measurements were performed us¬ 
ing the SExtractor package in dual image mode with the non- 
degraded K image adopted as the reference for source detection. 
The £lux_auto in all bands was measured with an aperture of 
2.5 times the Kron radius that includes > 96% of the total flux of 
the galaxy ( |Kron||T980 1 . Hereafter, the K-band magnitude was 
corrected to the total flux by measuring the growth curve of 
bright stars out to a radius of 8" (depending on the magnitude 
this correction ranges between 2%-4%). 

The space-based imaging from GALEX, IRAC and MIPS 
have more complicated PSP shapes and larger PWHM, there¬ 
fore photometry for these bands was performed separately (see 
Sec. 3.5 and 3.6 in |Muzzin et al.|2013| l. 

The photometry in all bands is corrected for Galactic dust at¬ 
tenuation using dust maps from |Schlegel et ak (1998 1 and using 
the Galactic Extinction Curve of |Cardelli et al.|p989[ ). The cor¬ 
rections were of the order of 15% in the GALEX bands, 5% in 
the optical and < 1% in the NIR and MIR. 

Star vs. galaxy separation was performed in the J - K versus 
u - J color space where there is a clear segregation between the 
two components (Pig. 3 of Muzzin et al.||2013] l. Sources were 
classified as galaxies if they match the following criteria; 


6. The galaxy sample 

To relate the occurrence of SN events to their parent stellar pop¬ 
ulation, we need to characterise the galaxy population in the sur¬ 
vey fields and in the redshift range explored by the SN search. To 
this aim, the extensive multi-wavelength coverage of both COS¬ 
MOS and CDPS gives a unique opportunity. In particular the 
analysis of deep multi-band surveys of the COSMOS field has al- 
ready been published (|Muzzin et al.|20131l, and we could retrieve 
the required information, such as photometric redshifts, galaxy 
masses and star formation rates directly from public catalogs. 
Por the CDPS fields we performed instead our own analysis but 

J - K > 0.18 X {u - J) - 0.75 for u - J < 3.0 

J - K > 0.08 x(u- J) - 0.45 otherwise (2) 

The photometric redshifts for the galaxv sample were ob- 

tained with the EAZT|''’| code (Brammer et al. 2008|l. EAZY 

fits the galaxy SEDs with a linear combinations of templates 
and includes optional flux- and redshift-based priors. In addi¬ 
tion, EAZY introduces a rest frame template error function to ac¬ 
count for wavelength dependent template mismatch. This func¬ 
tion gives different weights to different wavelength regions, and 
ensures that the formal redshift uncertainties are realistic. 

The template set adopted by|Muzzin et al. (2013 i includes: i) 

six templates derived from the PEGASE models (Fioc & Rocca- 

closely following the method described by Muzzin et al.|(2013|l. 

Volmerangel 199911, H) ^ rsd template from the models of Maras- 

In the following, we describe the detection and characterisation 

ton ( 2005]l, Hi) a 1 Gyr old single-burst (Bruzual & Chariot 2003|l 


of the galaxies in our search field. 


http://sngroup.oapd.inaf.it/asnc.html 


model to improve the fits for galaxies with post starburst-like 
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SEDs and in) a slightly dust-reddened young population to im¬ 


prove the hts for a population of UV-bright galaxies. Muzzin 
(2013|l chose to use the vl.O template error function, and 


et al. 


the K magnitude prior, and allowed photometric redshift solu¬ 
tions in the range 0 < z < 6. 

Photometric redshifts are extremely sensitive to errors in 
photometric zeropoints. A common procedure to address this 
problem is to rehne the zeropoints using a subsample of galax¬ 
ies with spectroscopic redshifts (e.g., Ilbert et al. 2006 [ Brammer 


|et al.|20rT] l. [Muzzin et al.| ( |2013]l used an iterative c ode devel¬ 
oped for the NMBS survey (see Whitaker et al.|20fT] l and found 
zeropoint offsets of the order of ~ 0.05 mag for the optical bands 
and of 0.1-0.2 mag for the NIR bands. 

To remain above the 90% completeness limit and guaran¬ 
tee the consistency with the CDFS catalog (see next section) we 
selected from the full COSMOS catalogue all galaxies with K 
band magnitude < 23.5. We further restrict the catalog to the sky 
area coverage of our held of view (1.15 deg^) and redshift range 
0 < z < 1 of interest for the SN search, obtaining a hnal count of 
67417 galaxies. 



Ks [ABmag] 


Fig. 8. Distribution of K magnitude for galaxies in CDFS (solid line) 
and COSMOS (dashed line) catalogues. 


6.2. VOICE-CDFS 

Areas of different sizes around the original CDFS field have been 
variously observed at different depths from the X-ray through 
the UV, Optical, IR to the Radio. The 0.5 deg^ Extended CDFS 
(ECDFS) multi-wavelength dataset has been carefully reduced 
and band-merged over the years (e.g. jCardamone et al.||2010[ 
[Hsu et aL]|2014| and reference therein). Conversely, most pub¬ 
lic multi-wavelength data over the VOICE-CDFS 4 deg^ area 
have been collected recently, and are not available as an homoge¬ 
neous database for our study. For our purposes we thus collected, 
merged and analysed most existing data ourselves. 

Available data over the VOICE-CDFS area include the fol¬ 
lowing: 

- GALEX UV deep imaging ( [Martin et al.|[2005b l. The 
GALEX photometry is from the GALEX GR6PIus7 data re- 
leasJ^ 

- SUDARE/VOICE u, g, r, i deep imaging (this work, Vaccari 
et al. in prep.) 

- VISTA Deep Extragalactic Observations (VIDEO, jJarvisj 
et al.|201^ Z,Y,J,H,K deep imaging 

- SERVS Spitzer Warm 3.6 and 4.5 micron deep imaging 
(jMauduit et al.|2012jl 

- SWIRE Spitzer IRAC and MIPS 7-band (3.6, 4.5 5.8, 8.0, 
24, 70, 160 micron) imaging ( [Lonsdale et al.|2003] l 

While data products are available as public catalog for 
most of the multi-wavelength surveys listed above, SERVS and 
SWIRE data were re-extracted and band-merged with all other 
datasets as part of the Spitzer Data Fusion project ( [Vaccari et al.| 
|2010[ http://www.mattiavaccari.net/df/). 

With the VIDEO survey still in progress, at the moment the 
sky areas covered by SUDARE and VIDEO do not fully over¬ 
lap: restricting our analysis to the overlapping region, for galaxy 
detection we lack a small portion of our VOICE-CDFS 1 and 
VOICE-CDFS2 (0.14 and 0.05 deg^, respectively) . However, 
for the estimate of SN rates in the cosmic volume surveyed by 
SUDARE we use the full area covered by the two fields. 

Deep image stacks have been obtained from SU¬ 
DARE/VOICE data as described in Sect. [2.1 j The VIDEO 

'* http://galex.stsci.edu/GR6/ 



redshifts 


Fig. 9. Redshift distribution (Zpeak) for galaxies in CDFS (solid line) and 
COSMOS (dashed line) catalogues. 


exposures were processed at the Cambridge Astronomical 
Survey Unit (CASU) using the pipeline developed specifically 
to the reduction of VIR CAM data, as py t of the VISTA Data 


Flow System (VDFS^ 

12 

Irwin et al. 

2004 

1 . The stacks produced 

by CASU were weightec 

mean com 

rined using SWarp. 


CDFS stacks in optical and NIR filters show a very small 
variation in seeing (ranging from 0.8"to 0.9") and we do not 
need to perform PSF homogenisation for measuring colours, but 
only to resample both VST and VISTA images to the same pixel 
scale of 0.21"pixeL' (for this we used SWarp). 

Source detection and photometry for VOICE-CDFS were 
performed with SExtractor in dual image mode with the K band 
image used as reference for source detection. 

The photometry was corrected for Galactic extinction corre¬ 
sponding to a flux correction of 3% in the optical and < 1% in 
the NIR. Then, we separated galaxies from stars using Eq. 


http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/data- 

processing 
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For all galaxies in the catalog we obtained photometric red- 
shift using the EAZY code adopting the same parameters and 
templates described in Sect. |6.1| The main difference between 
the two fields is the number of filters available for the analysis, 
12 filters for CDFS and 30 for COSMOS. To reduce catastrophic 
failures, we do not compute photometric redshifts for the sources 
that have detection in less than 6 filters (< 5%). 

Similar to Muzzin et al. ( 2013| l, the magnitude zero-points 
were verified using the iterative procedure developed by Bram- 
mer (p.c.). The procedure is based on the comparison of pho¬ 
tometric to spectroscopic redshifts: systematic deviations are 
translated into zero-point offsets corrections using K as the "an¬ 
chor" filter , the photometric scale is adjusted and EAZY re-run. 
We did not calculate offsets for GALEX and Spitzer bands. We 
found that g, r, i bands require small offsets (< 0.05 mag) while 
the u band requires an offset of 0.14 mag, and NIR bands of 
about 0.1 mag. 

Also for VOICE-CDFS we selected all galaxies with K band 
magnitude < 23.5 and redshift 0 < z < 1 that results in a final 
catalog of 92324 galaxies for VOICE-CDES that, allowing also 
for a small overlap of the two pointings, covers 2.05 deg^ . 

The distribution of K band magnitudes and photometric red¬ 
shifts for the COSMOS and VOICE-CDFS galaxy samples are 
shown in Fig. and in Fig. [9]respectively. 


6.3. Accuracy of photometric redshifts 

We explored different methods to assess the quality of photomet¬ 
ric redshifts : /) analysing the width of confidence intervals, and 
quality measurements provided by EAZY, ii) comparing differ¬ 
ent redshift estimators. Hi) comparing the photometric redshifts 
with available spectroscopic redshifts, iv) comparing our esti¬ 
mates with photometric redshifts from other groups. 


6.3.1. Internal error estimates 

EAZY provides multiple estimators of the photometric redshifts 
among which we choose Zpeak corresponding to the peak of the 
redshift probability distribution F(z). As a measure of the uncer¬ 
tainty, the code provides 68, 95 and 99% confidence intervals 
calculated by integrating the P{z). The confidence intervals are a 
strong function of the galaxy apparent magnitude and redshift, as 
shown in Eigjl^for the 68% level. The narrower confidence in¬ 
tervals for the COSMOS field with respect to the VOICE-CDFS 
field are due to the better sampling of the SED for the galaxies 
of the former field. 

EAZY provides also a redshift quality parameter, eO that 
is intended as a robust estimate of the reliability of the photo¬ 
metric redshift ( [Brammer et al.|[2008 1 . Poor fits (Q^ > 1) may 
be caused by uncertainties in the photometry, poor match of the 
intrinsic SED from the adopted templates, or degeneracies and 
nonlinear mapping in the colour-z space. We found good qual¬ 
ity photometric redshifts (Q^ < 1) for 75% and for 93% of the 
galaxies in CDFS and COSMOS field, respectively. 

In several cases the P(z) function is multimodal, so that Zpeak, 
that corresponds to the peak of P{z), does not properly reflect 
the probability distribution. |Wittman| ( |2009| l introduced a very 
simple alternative estimator that represents the redshift proba¬ 
bility distribution, incorporating the redshift uncertainties. This 


Qz - Na-i p” f 2 number of photometric mea¬ 
surements used in the fit, is the 99% confidence intervals and 

Paz=o.2 is the fractional probability that the redshift lies within ±0.2 of 
the nominal value. 
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Fig. 10. The width of the 68% confidence intervals computed from the 
redshift probability distribution as a function of galaxy magnitude (top 
panel) and redshift (bottom panel). 


redshift estimator is drawn randomly from the P(z) and denoted 
with zmc because it results from a Monte Carlo sampling of the 
full P(z). The difference between Zpeak and zmc can be used as an 
indication of the internal uncertainties of photometric redshifts. 
The difference in the redshift distribution obtained with different 
redshift estimators can be seen in Fig. [m 


6.3.2. Comparison with Spectroscopic Redshifts 

Spectroscopic redshifts are available for a fairly large number 
of galaxies for both our fields. The spectroscopic redshifts for 
4733 galaxies in the COSMOS field were taken from Muzzin 
|et al.| ( |2013| l while for the CDFS field the data for 3362 galax¬ 
ies were collected from the literature, from different sources and 
with different quality flags. A comparison of photometric and 
spectroscopic redshifts for these subsample is shown in Fig. fl^ 
We calculated the normalized, median absolute deviatiotyg 
(NMAD) which is less sensitive to outliers compared to the 


Brammer et al. 


2008 


o-nmad = 1-48 X median as in 

I ^+Zspec I 

where Az = (Zpimt — Zspec). The normalization factor ot i.48 ensures that 
NMAD of a Gaussian distribution is equal to its standard deviation, and 
the subtraction of median)Az) removes possible systematic offsets 
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Fig. 11. The distribution of Zpeak (solid line) and zmc (dashed line) of 
the galaxies in COSMOS (top panel) and in CDFS catalogues (bottom 
panel). 


standard deviation ( [Brammer et al.|2008| l. For CDFS we found 
c^NMAD = 0.02 which is comparable to that of other surveys 
with a similar number of filters, whereas for COSMOS (Tnmad = 
0.005. 

Another useful indication of the photometric redshift quality 
is the fraction of "catastrophic" redshifts defined as the fraction 
of galaxies for which |zphot - Zspec| /(I +Zspec) > Sctnmad- For the 
CDFS field we found a fairly large fraction of catastrophic red¬ 
shifts (~ 14%). After removing these outliers, the rms dispersion 
Az/(1 +z) = 0.02. The same analysis for the COSMOS field (see 
Muzzin et al.||2013] for details) gives a fraction of 5cr outliers 
as low as 4% and a very small rms dispersion for the rest of the 
sample (0.005). 

6.3.3. Comparison with Zphot from other surveys 

The comparison between photometric and spectroscopic red¬ 
shifts is biased towards brighter galaxies for which it is easier 
to observe the spectrum. To analyse the accuracy of our photo¬ 
metric redshifts in a wider luminosity range, we compare our 
estimates to those obtained by the Multiwavelength Survey by 
Yale-Chile (MUSYC, [Cardamone et al.|2010] l which covers the 
~ 30' X 30' “Extended” Chandra Deep Field-South (that is in¬ 
cluded in CDFSl) with 18 medium-band filter optical imag¬ 
ing from the Subaru telescope, 10 broad-band optical and NIR 
imaging from the ESO MPG 2.2 m (Garching-Bonn Deep Sur¬ 
vey), ESO NTT and the CTIO Blanco telescopes along with 4 
MIR bands IRAC imaging from Spitzer SIMPLE project. The 
MUSYC catalog lists BVR-selected sources with photometric 
redshifts derived with the EAZY program. Therefore the main 
difference is that the MUSYC catalogue makes use of a much 
larger number of filters compared with SUDARE, which signifi¬ 
cantly improves the photometric redshift accuracy. 

Cross correlating the two catalogs with a search radius of 2" 
we found 1830 common galaxies. In Eigure[^ we plot the dif¬ 
ferences between the Zphot estimates as a function of the 
We find evidence of a some systematic difference at low red¬ 
shifts z < 0.3, with the Zphot from SUDARE being higher, but 
overall the two catalogs shows a fair agreement with a scatter 
Az/(1 +z) - 0.05 and a 5cr outlier fraction of 10%. 
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Fig. 12. Comparison of photometric vs. spectroscopic redshifts from 
COSMOS (top panel) and CDFS (bottom panel). 


7. Computing SN rates 

To compute the SN rate we need to introduce the method of the 
control time (CT, |Zwicky|[l942[ ). The CT of one observation is 
defined as the interval of time during which a SN occuri'ing at 
a given redshift is expected to remain above the detection limit 
of the image. The total CT of an observing campaign is prop¬ 
erly computed by summing the CT of the individual observa¬ 
tions ( |Cappellaro et al.|1997 1 . Then, the SN rate is computed as 
the number of events detected in the survey divided by the total 
CT. 

The CT depends on the SN luminosity and light curve evo¬ 
lution and is therefore different for different SN types. We con- 
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Fig. 13. Comparison of our photometric redshift for CDFS and MUS YC 
photometric redshifts. 


sidered separately the following main SN types: la, Ib/c, II (in¬ 
cluding IIP and IIL), Iln and SLSN. 


7. 1. The control time 

To compute the CT we select a template light curve represen¬ 
tative of a given SN subtype (SNi), a redshift (z, in the range 
0 < z < 1), and an extinction value in the range Ay = 0 - 2 mag 
(in the host galaxy rest frame). To take into account the diversity 
of the photometric evolution for SNe of different types we used 
a wide collection of light curve templates (listed in Tab. [^. We 
considered four representative subtypes for thermonuclear SNe 
(normal, bright, faint and peculiar) , six subtypes for hydrogen 
rich SNe (IIP, IIP faint, IIL, Ilb, Iln, plus peculiars) three sub- 
types for stripped envelope SNe (Ib, Ic, Ic broad line) along with 
a template for SLSN. In some cases, we use a few templates for 
the same SN subtype to take into account the photometric vari¬ 
ance within the class. 

Then: 


- we define a useful range for the epochs of explosion. In order 
to be detectable in our search, a SN needs to explode in the 
interval [fo - 365^f, f/^], where to and tj^ are respectively the 
epochs of the first and last observations of the given field. 
In fact, for the redshift range of interest of our survey, a SN 
exploded 1 yr earlier than the first observation is far too faint 
to be detected. 

- we compute the expected magnitude, m,, at each epoch of 
observations, f, (with i = 1,2, ....K, where K is the number 
of observations), for a SN that explodes at an epoch xj in¬ 
cluded in the time interval defined above. To derive these 
estimates, we use the SN template light curve, the proper 
K-corrections, the distance modulus for the selected redshift 
and the adopted extinction. 

- the detection probability pi{xj) of the simulated event at each 
observing epoch is given by the detection efficiency for the 
expected magnitude, e,(OT/), estimated as described in Sect[^ 
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The detection probability for the whole observing campaign 
is derived as the complement of the probability of not detec¬ 
tion at any of the epochs, that is p{,Xj) - 1 -R^o R “ 

- we simulate a number N of events exploring the possible 
epochs of explosion, in the interval [fo - 365d, t^]. We can 
then compute: 


CTsNi.Esvfe) - (^K - to + 365) 


N 


(3) 


where f is expressed in Julian Day. The accuracy of the CT 
computation above depends on the sampling for the explo¬ 
sion epoch in the defined interval. After some experiments, 
we found that a sampling of Id is more than adequate, con¬ 
sidering also the contribution of other error sources. 


- for the extinction distribution, following Neill et al. ( 2006) 1 
we adopted an half-normal distribution with (Te(b-v) - 0-2. 
We adopt the same distribution for all SN subtypes although 
we may expect that different SN types, exploding in differ¬ 
ent environments, may suffer different amount of extinction. 
In particular, the distribution of Neill et al. (2006 1 was de¬ 
rived for SN la and is likely to underestimate the effect for 
CC SNe. In Sect. |7.3| we verify (a posteriori) the consistency 
of our assumptions about the extinction distribution and esti¬ 
mate how its uncertainty propagates in the systematic uncer¬ 
tainty of SN rates. 


- finally, the CT for each of the main SN types was com¬ 
puted by accounting for the subtype distribution and for the 
adopted extinction distribution (details below): 


CTsn(z) = zz /sNi SEbv CT SNi,EBv (^) ■ (4) 

SNi Egv 

where fsNi is the SN subtype fraction and ggBv the distribu¬ 
tion of colour excess E(B-V) (multiple templates for a given 
subtype are given equal weight). 


In principle, with sufficient statistics and accurate subtype 
classification, the fractional contribution of the different sub- 
types can be derived from the distribution of detected events. 
However, at the current stage of the project, the event statistics is 
not large enough and we adopted the subtype distribution from 
the literature. In particular, we adopted the fractions of different 
SN subtypes obtained by Li et al. (201 lb i with the exception 


of the fraction of faint type II SNe that is from Pastorello et al. 


(2004 1 while for bright Ic SNe we refer to Podsiadlowski et al. 
(2004 1 . Our adopted subtype distribution is: 


- type la: 70% normal, 10% bright 1991T-like, 15% faint 
1991bg-like and 5% 2002cx-like ; 

- type II: 60% IIP, 10% 2005cs-like , 10% 1987A-like, 10% 
IIL and 10% Ilb, 

- type Ib/c: 27% Ib, 68% Ic and 5% 1998bw-like, 

- type Iln: 45% 1998S-like, 45% 2010jl-like, 10% 2005gj-like, 


We stress that these subtype distributions are obtained from a 
local sample and it is possible, or even expected, that they evolve 
with cosmic time. In Sect. |7.3| we estimate the uncertainty im¬ 
plied by this assumption. 
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Fig. 14. Control time as a function of redshift for different SN types (in 
days per deg^) averaged across the three survey fields. 


7.2. SN rate per unit volume 

The volumetric SN rates per redshift bins in the range 0 < z < 1 
is calculated as: 


_ (1 +z) Nsn(z) 
Viz) crsNfe)' 


(5) 


where Nsn(z) is the number of SNe of the given type in the spe¬ 
cific redshift bin, C7’sn(z) is the control time and the factor (1 -Hz) 
corrects for time dilation. Viz) is the comoving volume for the 
given redshift bin that is computed as: 


V(z) = 


Ttt 0 


n dz' 

3 41253 

Ho J 

zi sjQ.n^i l + z')^ . 


Mpc^ 


(6) 


where 0 is the search area in deg^ and z is the mid-point of the 
redshift bin with extremes zi, Zi- 


7.3. Statistical and systematic errors 

We derive the 1-cr lower and upper confidence limits from the 
event statistics as in |Gehrels] ( |1986| l. Then, these values are con¬ 
verted into confidence limits of SN rates through error propaga¬ 
tion of Eq.|^ 

There are many sources of systematic errors. To estimate 
each specific contribution we performed a number of experi¬ 
ments calculating SN rates under different assumptions. 

Transient Misclassification 

As described in Sect |4.H for a fraction of SN candidates 
(23%) the SN confirmation remains uncertain. These PSNe 
are attributed a weight 0.5 in the rate calculation. To obtain an 
estimate of the impact of this assumption, we compute the rate 
in the extreme cases assuming for these events a weight 0 and 
1 respectively. As an error estimate, we take the deviation from 
the reference value of the rates obtained in the two extreme 
cases. It turns out the error is of the order of 10 - 15% (Tab.|^. 
One concern is that we set arbitrary thresholds for P ^2 and Apt to 
attribute the flag of PSN. To test the impact of this assumption, 
we computed the SN rate adopting different thresholds, namely 


10 ^ or 10 ® for P^i and 5 or 9 for Apt. We found that in all 
cases the deviations for the reference value are < 10% (typically 
~ 5%). 


SN photometric typing 

For the errors of on SN typing we adopt the values discussed 
in Sect. I^that is 10% for type la, 25% for type II and 40% for 
type Ib/c and Iln, independently on redshift. It appears that the 
error on SN typing has in general a moderate impact for type 
la, and instead it is one of the dominant sources for SN CC in 
general. 

Subtype distribution 

The adopted SN subtype distribution affects the estimate of 
SN rates because the subtypes have different light curves hence 
different control times. We consider an error of 50% for the 
fraction of the subclasses and, as an estimate of the contribution 
to the systematic error, we take the range of values of the SN 
rates obtained with the extreme subtype distribution. This is a 
significant source of error, typically 10-20% but with a peak of 
40% for type Iln SNe. 

Detection efficiencies 

We performed MonteCarlo simulations assuming that the 
value of the detection magnitude limit of each observation has 
normal error distribution with cr - 1.0 mag. We found that 
the frequent monitoring of our survey means that the large 
uncertainty in the detection efficiency for each single epoch 
does not have a strong impact on the overall uncertainty. The 
propagated error on the rates is < 10%. 

Plost galaxy extinction 

In our computation we adopt a half-normal distribution of 
EiB - V) with cr = 0.2 mag for both type la and CC SNe. 
To estimate the effect of this assumption, the SN rates have 
been recalculated assuming a distribution with cr = 0.1 mag 
and cr = 0.3 mag. We evaluate that the error on the rates is of 
the order of 5 - 10%. The uncertainty is more critical for CC 
SNe (11%) and for the highest redshift bin of type la SNe (14%). 

The consistency of the adopted extinction distribution was 
verified a-posteriori. We computed estimates of the SN rates 
for a range of crE{B-v) values ranging from 0 to 0.5 mag. For 
each adopted (Te(b-v), we computed the expected distribution 
of extinction of the detected SNe. This is different from the 
intrinsic distribution because of the bias against the detection of 
SNe with high extinction which have a shorter control time. The 
expected extinction distribution is compared with the observed 
distribution (Fig. [T5] l and the best matching cr is determined 
using a Kolmogorov-Smirnov 2-side test. We found a best 
match for (Te{b-v) - 0.25 mag using the full SN sample or 
crE(B-v) - 0.28 mag using only the type la events and excluding 
probable SNe (PSN). Given the uncertainties, we consider 
that this values are consistent with the adopted distribution 
from Neill et aT| ( |2006| l. We remind that the adopted EiB - V) 
distribution was used only for the CT calculation and not for the 
SN photometric classifier. 


Photometric redshifts 

We compare the results obtained using the two alternative 
photo-z estimatator Zpeak and Zmc (c.f. Sect 6.3. li. It turns out 


that this is the most significant source of error, especially for 
type la SNe. A detailed analysis shows that the most important 
effect is for the redshift of SN host galaxies, while the effect on 
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Fig. 15. Predicted (line) vs. observed (bars) cumulative extinction dis¬ 
tribution. 


Table 6. Relative systematic errors 


< z > 

0.10 

la 

0.25 0.45 

0.65 

CC 

0.10 0.25 

lln 

0.55 

PSN 

- 

0.04 

0.11 

0.12 

0.26 

0.16 

0.23 

SN tvping 

0.10 

0.10 

0.10 

0.10 

0.32 

0.31 

0.40 

subtype distr. 

0.08 

0.10 

0.20 

0.24 

0.07 

0.12 

0.40 

detection eff. 

0.02 

0.05 

0.05 

0.07 

0.02 

0.05 

0.03 

extinction 

0.04 

0.03 

0.07 

0.14 

0.02 

0.11 

0.05 

z distribution 

0.33 

0.21 

0.17 

0.21 

0.09 

0.07 

0.03 

all 

0.36 

0.26 

0.31 

0.39 

0.43 

0.39 

0.61 


the control time of the galaxy population has a smaller impact. 

Cosmic variance 

The possible under/over-density of galaxies in the held of 
view due to cosmic variance impacts on the SN rate measure¬ 
ments. Using the cosmic variance calculator of 
( |2008| l we found that cosmic variance can add 
the SN rate of 5-10% but for the low redshift bin the variance 
can be as large as 15-20%. We note that the cosmic variance 
bias is averaged out when rate measurements from different sky 
helds are analyzed together as for SUDARE. 

In Table we report the individual systematic errors along 
with the overall error obtained by their sum in quadrature. We 
do not include the effect of cosmic variance in the etTor budget, 
since this is not a measurement error. Rather, this is an uncer¬ 
tainty related to the particular galaxy sampling in our survey. 

The overall systematic etTor is typically of the order of 30 - 
40%, and is larger than the statistical error. 


Trent! & Stiavelli 
an uncertainty of 


8. SN rates as a function of cosmic time 

Our SN rates per unit volume are reported in Tab [7] Columns 1 
and 2 report the SN type and redshift bin, col. 3 gives the num¬ 
ber of SNe (the number of PSNe is in parenthesis), col. 4 the rate 
measurements, cols. 5 and 6 the statistical and systematic errors 
respectively. The redshift bins were chosen to include a signif¬ 
icant number of SNe (a minimum number of 10 SNe) with the 
exception of the nearest redshift bin (0.05 < z < 0.15) , where 
we collected only few SNe. 
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Table 7. SN rates per unit volume [10^“* yr“* Mpc“^]. Upper limits were 
computed for a reference value of 3 events. In this case, from Poisson 
statistics, the probability of obtaining a null results is < 5%. 


SN type 

Zbin 

SNe 

rate 

stat. 

syst. 


0.05-0.15 

3.0 

0.55 

-0.29 -hO.50 

+0.20 

la 

0.15-0.35 

12.7 

0.39 

-0.12 -hO.13 

+0.10 

0.35-0.55 

23.0 

0.52 

-0.13 -hO.II 

+0.16 


0.55-0.75 

17.4 

0.69 

-0.18 -hO.19 

+0.27 


CC 

0.05-0.15 

5.9 

1.13 

-0.53 +0.62 

+0.49 

0.15-0.35 

26.2 

1.21 

-0.27 +0.27 

+0.47 

II 

0.15-0.35 

13.4 

0.69 

-0.18 +0.16 

+0.24 

Ib/c 

0.15-0.35 

9.3 

0.48 

-0.17 +0.19 

+0.23 

Iln 

0.15-0.35 

3.5 

0.043 

-0.026 +0.030 

+0.026 

0.35-0.75 

5.8 

0.017 

-0.009 +0.009 

+0.010 

SLSN 

0.35-0.75 

- 

< 0.009 




8.1. Core collapse SNe 


Fig. 16 shows a comparison of our estimate of the rate of CC 


SNe with all measurements available in the literature. To obtain 
the CC SN rate we cumulated type 11, Ib/c and lln events. 

Our results are in good agreement with other measurements. 
Note that in Fig. T^we re port the value of Melinder et al. (2012i 
and pahlen et al.|(|2012ll w ith no correction for the fraction of 
hidden SNe ( Mattila et al.||201^ . We will return on this point 
later. 

Given the short lifetime of their progenitors (< 30 Myr), 
there is a simple, direct relation between the CC SN and the cur¬ 
rent SF rate: 


rcciz) = Kcc X ij/(z) (7) 

where il/(z) is the SFR and Kcc is the number of stars per unit 
mass that produce CC SNe, that is: 


Kcc 


rmu.cc 
Jmccc ^ 

m<f)(m)dm 


( 8 ) 


where (f>(m) is the initial mass function (IMF), me and mu are 
the extreme limits of the stellar mass range, mL.cc and »*u,cc the 
mass range of CC SN progenitors. 

Assuming that Kcc does not evolve significantly in the red¬ 
shift range of interest, the evolution of the CC SN rates with 
redshift is a direct tracer of the cosmic SF history (SFH). Con¬ 
versely, we can use existing estimates of the SFH to compute 
the expected CC SN rate, assuming a mass range for their pro¬ 
genitors. In order to do this consistently one has to use the same 
IMF (or Kcc) adopted to derive the SFR. Indeed, although in 
Eqj^Kcc depends on the IMF, the ratio between the cosmic SFR 
and CC rate does not give a real indication on the IMF, since 
both quantities actually trace the number of massive stars, which 
produce both UV photons and CC SN events. The formal depen¬ 
dence on the IMF of this ratio is introduced by the extrapolation 
factor used when deriving the SFH from luminosity measure¬ 
ment to convert the number of massive stars formed at the vari¬ 
ous redshifts into the total stellar mass formed. 
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Fig. 16. CC SN rate per unit volume. All measurements do not account for the correction for hidden SNe. To obtain the predicted SN rate from 
the measured SFR we adopt 8, 40 Mq as the lower and upper mass limits for SN CC progenitors and the proper IMF, Salpeter for |Madau &| 
|DickinsoT|j2014[l and SalA for |Hopkins & BeacomH2006) . The dashed lines show the predicted SN rate assuming the fraction of hidden SNe 
given in (|Mattila et al.|2012|>. 


The CC progenitor mass range is still uncertain, both for the 
low and upper limit. Stellar evolution models suggest a typical 
range of 9 - 40 Mq ( |Heger et al.|2003| l for CC SNe, though the 
upper limit strongly depends on metallicity and other factors, 
e.g. rotation or binarity. In recent years, for a number of nearby 
CC SN it was feasible to search for the progenitor star in archival 
pre-explosion images( [Smartt|20()9||2015[ and reference therein). 
This allows to obtain an estimate of the masses of their progen¬ 
itor stars or, if not detected, of upper limits. By comparing the 
observed mass distribution with the IMF, it was argued that the 
minimum initial mass is 8 + 1 Mq. The same analysis also sug¬ 
gests a paucity of progenitors of SN II with mass greater than 
20 Mg, which would indicate that these stars collapse directly to 
a black hole, without producing a bright optical transient ( [Smartt 
|2009| l. However this result needs to be confirmed and hereafter, 
following the trend of the literature in the field, we adopt an up¬ 
per limit of 40 Mg. 

With a mass range 8 - 40 Mg for the SN CC progenitors 
we obtain a scale factor Kqc - 6.7 x 10“^ Mg^ for a standard 
Salpeter IMF or Kqc = 8.8 x 10~^ Mg' for a modified Salpeter 
IMF (SalA), with a slope of -1.3 below 0.5 Mg (similar to what 
adopted in|Hopkins & Beacom|2006|. 


It has been claimed that assuming the 8-40 Mg mass range. 


the comparisons between the SFH from Hopkins & Beacom 
P006| l (hereafter HB06) and the published measurements of 
CC SN rates showed a discrepancy of a factor 2 at all red- 
shifts ( |Botticella et al.|2008[ |Bazin et al.|2009| l. |Horiuchi et'af] 
( |2011| l argued that this indicated a "supernova rate problem" for 
which they proposed some possible explanations: either many 
CC SNe are missed in the optical searches because of heavy 
dust-obscuration or there is a significant fraction of intrinsically 
very faint (or dark) SNe when, after core collapse, the whole 
ejecta falls back onto the black hole. 


On the other hand, |Botticella et al.| ( |2012T l found that CC SN 
rate in a sample of galaxies within 11 Mpc is consistent with that 


expected from the SFR derived from FUV luminosities. Taylor 


|et al.| ( [20T4| l based on the SDSS-II SN sample estimated that the 
fraction of missing events is likely of the order of 20%. |Gerke| 


et al. (2014 1 performed a search for failed SNe monitoring a sam¬ 


ple of nearby galaxies (< 10 Mpc). After 4 yr they found only 
one candidate suggesting an upper limit of 40% for the fraction 
of dark events among CC SNe that, unfortunately is not yet a 
strong constraint. 
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To detect the CC SNe hidden by strong extinction, a few 
infrared SN searches have been performed in local starburst 
galaxies (Maiolino et al. 20021 Mannucci et al. 2003[ [Mattila 
|& Meikle||2001[ Miluzio et al.(|2013|l in some cases exploiting 


adaptive optics (|Cresci et al.|2()07[ Mattila et al.|2007] Kankare 


et al.| 2008 2012| l to improve the spatial resolution. However, 
despite the efforts, it was not possible to unveil the hidden SNe. 

An alternative approach to estimate the fraction of hidden 
SNe was made with the conservative assumption that all SNe 
in the nucleus of luminous and ultra-luminous infrared galaxies 


(LIRGs and ULIRGs) are lost by optical survey ( Mannucci et al.| 
2007| l. With this approach, Mattila et al. ( |2012 suggested that 
the fraction of missed SNe increases from the average local value 
of - 19%to~ 38%atz~ 1.2. 

The CC SN rate predicted using two different SFH from 
HB06 and the recent results of |Madau & Dickinson| ( |2014l l (here¬ 
after MD14) are shown in Fig. 16 The two SFH lead to differ¬ 
ent predictions with a discrepancy of about a factor 2. Indeed, 
the cosmic SFR derived by MD14 is lower than the HB06’s vir¬ 
tually at all redshifts. In addition, the MD14 SF rates assume a 
straight Salpeter IMF, so that the number of massive stars formed 
in the Universe at all epochs is further diminished, when com¬ 
pared to predictions obtained with HB06’s SFH. Both factors 
concur in producing the final result shown in Fig{^ The predic¬ 
tions based on the MD14 SFH are in good agreement with the 
data. We note that by applying the Mattila et al.| correction for 


hidden SNe (dashed green lines in the figure) improves the fit at 
low redshift but gives a worse comparison at high redshift. 

The HB06 SFH instead over-predicts the CC rate, if the pro¬ 
genitors come from the mass range 8 to 40 Mq. In this case, cor¬ 
recting for hidden SN improves the agreement at high redshift, 
but still overestimates the CC SN rates at z < 0.4 (dashed blue 
line). 

On the other hand, as we mentioned above, the uncertainty 
in the CC progenitor mass range is still significant. Indeed, most 
recent data collection seems too indicate a lower mass limit as 
high as 9 - IOMq ( Smartt|2015 i. If we use this in combination 
with the upper mass limits for CC SN progenitors of ~ 20 Mq 
( e.g. |Smartt|[2009[ |Gerke et al.||2014| and also include the cor¬ 
rection for hidden SNe, it would result that the SFR severely 
under-predict the CC rate (by over a factor 3 if we refer to the 
MD14SFR). 

All together, it is fair to say that the both statistic and system- 
atics errors on SN rate and SFR measurements and the uncertain¬ 
ties on the progenitor mass range are too large to invoke a "SN 
rate problem" and hence to speculate on possible explanations. 

One of the goals we aim to achieve with our survey it to 
obtain measurements of the evolution of specific SN subtype. 
While the statistics of the present sample is still small, we can 
however obtain some preliminary measurements. 


8.2. SN lln 

Our estimate of the rate of type Iln SNe is uncertain for two 
reasons: SNe Iln are rare and the event statistics is very poor. 
In addition, the variety in luminosity and light curve evolution 
that, in some cases, mimic those of other SNe (e.g. SN IIL or 
SLSN), makes the photometric classification very uncertain (cf. 
Sect. |7.3| l. However, due to the intrinsically bright and slowly 
evolving luminosity we could detect type Iln SNe in a redshfit 
range comparable to that of SN la. 

In the 0.15-0.35 redshift bin, we estimated that type Iln are 
4 + 3% of all CC SNe. This number is consistent with the 6 + 2% 
value measured in the local Universe (|M etal.|2011a| l. 

On the other hand, the apparent decrease of the rate at higher 
redshift (a factor ~ 2.5 in the redshift bin 0.35-0.75 compared 
with the nearest bin) appears at odd considering that the overall 
CC SN rate in the same redshift interval increases by about the 
same factor. 

Either there is a strong evolution of the type Iln rate with 
redshift or, in our search, we are missing (or mis-classifying) 
over 2/3 of the distant type Iln. Both explanations are difficult to 
accept: we will need to verify this result at the end of our survey 
with better statistics and, possibly, improved template list. 

8.3. Superluminous SNe 

Super-luminous supernovae (SLSNe) that radiate more than 
10 ^"^ergs“' at their peak luminosity (about 100 times more lu¬ 
minous than usual Type la and CC SNe) have recently been dis¬ 
covered in faint galaxies typically at high redshifts. The origin 
of these events is still unclear: the host environment and energet¬ 
ics suggest massive stellar explosions but their power source is 
still a matter of debate. In fact, several sub-classes have been in¬ 
troduced possibly related to different explosions scenarios (|Gal-| 
|Yam|2012| l. 

We did not detect any SLSN in our surveyed volume up to 
z = 0.75. However the null result for SLSN can provide in¬ 
teresting constraints on the rates of these objects. Lrom sim¬ 
ple Poisson statistics we find that the probability of obtaining 
a null result is 5% when the expected values is 3.0. Therefore 
one should expect that the rate of SLSNe is no higher than 
9 X lO^^yr 'Mpc“^ at a mean redshift z ~ 0.5. This firm up- 


ner limit is consistent with the rate estimated bv 

Quimby et al. 

(j2013 

that IS 

1 of 2.0 ^q g X 10 ^ yr 'Mpc ^ at a mean rec 
1 SLSN for each 500 CC SNe. 

shift z = 0.16 


We note that between z = 0.15 and z = 0.5 the CC rate 
measurements increases by almost a factor 3. Our upper limit 
does not preclude a similar increase for the SLSN rate; clearly 
we need to obtain a more significative result that will become 
feasible when our survey is completed. 


8.1.1. SN Ib/c 


8.4. SNe la 


We found that at the mean redshift z = 0.25, type Ib/c are 40 + 
13% of CC SNe. This compares very well with the estimates 
of Li et al. ( 201 la[ ) for the local Universe, measuring a fraction 
of Ib/c that ranges from 46 + 17% in early spiral galaxies to 
20 + 5% in late spirals, with an average value of 33 + 9%. The 
physical reason of the difference in CC SN population in early 
and late spirals is not understood, though it is possibly related to 
a metallicity effect ( |Li et al.|2011a[ ). 

Given the limited statistics of our sample, we can only con¬ 
clude that there is no evidence for evolution with redshift of the 
Ib/c fraction. 


Our measurements of the rate of SN la are shown in Lig. [TT] 
along with all those available from the literatur^^ 

Our results appear in agreement with other measurements 
within the statistical errors. It may be noticed, however, that our 


We do not plot original measurements that have been later revised or 


superseded as follow :|Madgwick et al.|(|2003) by|Graur & Maoz 

201311, 

Poznanski et al.H2007b|l by Graur et al. 

([201 l[l. 

Dahlen et al.[ 

2004|; 

Kuznetsova et al.|(|2008J by 

Dahlen et al. 

12008J, I 

Barris & Tonry[2006J 

by [Rodney & Tonry|p010 

1 , [Neill et al. 

([2006[[2007J bylPerrett et al. 

i2012|l 
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Table 8. Average of SN la rate measurements per redshift bin (units 
10-“* yr-‘ Mpc^^). 


^bin 

< z > 

ria 

(T 

N 

0.00-0.15 

0.05 

0.25 

0.05 

6 

0.15-0.35 

0.25 

0.29 

0.07 

7 

0.35-0.55 

0.45 

0.44 

0.11 

9 

0.55-0.75 

0.65 

0.58 

0.14 

8 

0.75-1.00 

0.84 

0.64 

0.20 

11 

1.00-1.50 

1.16 

0.87 

0.22 

7 

1.50-2.00 

1.64 

0.63 

0.22 

5 


4 

2 


4 

2 


4 

2 


0 

Fig. 18. Histogram of published estimates of the SN la rate in the three 
redshift bins of our measurements. All values are scaled to the mean 
redshift of the bin assuming a linear evolution of the rate with redshift 
(see text). In red we show our measurements. The black lines are the 
Gaussian curve whose mean and variance are computed from the data 
and reported in each panel’s legend. The averages and dispersions were 
computed by weighting the individual measurements with the inverse 
of their statistical errors. 



<z > =0.65 
<r-i^> =0.58 ±0.15 



SN rate 


estimates seem on the high side compared with the bulk of pub¬ 
lished measurements. 

In Fig. we plot the histogram of the measurements for the 
three bins of redshifts corresponding to our measurements. The 
effect of the rate evolution within each of these redshift bins has 
been removed by scaling the measurements to the mean redshift 
of the bin assuming that the rates scales as rja oc 0.6 x z, that, as 
a first order approximation, fits the rate evolution up to redshift 
~ 1 (the exact slope of the relation may be slightly different but 
it is not crucial for the comparison we are performing here). 

It appears from the Figure that for each bin the distribution of 
measurements (our own included) are consistent with a normal 
distribution and the dispersion is well understood considering 
the statistical and systematic errors affecting the measurements. 

Given that, in the following we will use average values as the 
best estimates of the SN la rate for the comparison with models. 
The average rates per redshift bin are reported in Tab. where 
col. 1 gives the redshift bin, col. 2 the average redshift, col. 3 
and 4 the average rate and dispersion and col. 5 the number of 
measurements per bin. Notice that for redshift z > 0.75 in the 
computation of the average rate we did not correct the individual 
measurements for the possible rate evolution inside the bin. 



Fig. 19. Distribution functions of the delay times selected for our the¬ 
oretical predictions for the Single Degenerate (green), and Double De¬ 
generate models (blue and red, see text for more details). The thin lines 
show the cumulative fraction of events as a function of time 


In Fig. 1^ we compare the average SN la rate measurements 
with the expected evolution for different progenitor scenarios 
predicted by Greggio (20051. 

Models in |Greggiof ( 2005 1 ) assume that SNIa progenitors are 
close binary systems which attain explosion upon reaching the 
Chandrasekhar mass either because of mass accretion from a 
companion star (single degenerate, SD) or of merging with an¬ 
other WD (double degenerate, DD). The delay between the birth 
of the binary system and its final explosion ranges from ~ 40 
Myr to the Hubble time, so that at each epoch the SN events in a 
galaxy are the result of the contributions of all past stellar gener¬ 
ations. Following Greggio (20051, the expected SN la rate at the 
time t is: 


X min(t,Tx) 

/la(T)l/^(f - T)dT (9) 

where Kia is the number of SN la progenitors per unit mass 
of the stellar generation, /ia(T) is the distribution function of the 
delay times and il/{t - t) is the star formation rate at the epoch 
t-T. The integration is extended over the full range of the delay 
time T in the range r, and min(t, t^), with t, and being the min¬ 
imum and maximum possible delay times for a given progenitor 
scenario. According to stellar evolution, f\Jj) is a decreasing 
function of the delay time, with a slope which depends on de¬ 
tails of the scenario leading to the SN explosion. One can then 
use Eq.|^to constrain the progenitor’s model using the trend of 
the SNIa rate with cosmic time, after specifying the cosmic SFH. 
In the following we adopt the MD14 SFH and assume that Ki^ 
does not vary with cosmic time. 

We select three DTD models, plotted in Fig.[T^ and test their 
predictions for the cosmic SNIa rate. The models include a sin¬ 
gle degenerate realisation (SD), and two flavours of the double 
degenerates, either with a close (DDC) or wide (DDW) binary 
separation predicting a steep and a mildly decreasing distribu¬ 
tion of the delay times, respectively (see|Greggio|2005||2010| for 
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Fig. 17. Our estimates of the SN la rate at z = 0.25,0.45,0.65 are compared with the other values from literature. The rate of |Cappellaro et ah] 
( |1999^ ; [Hardin et ^ ( |2000| l; |Madgwick et al.|p003| l; [Blanc et al.| ( |2004^ were given per unit luminosity. They are conv erted in rate per unit vo lume 


using the following relation of the luminosity density as a function of redshift: 7 b(z) = (1.03 + 1.76z) x 10° L® Mpc ° (Botticella et al. 


measurements of [Perrett et al.|(|2012|> are scaled up by 15% to account for the fact that they not include the faint SN 1991bg-like events. 


2008 1 . The 


more details). The selected models correspond to a very different 
time evolution following a burst of star formation. For the DDC, 
SD and DDW models, 50% of the explosions occur within the 
first 0.45, 1 and 1.6 Gyr respectively, while the fraction of events 
within 500 Myr is 0.55, 0.3 and 0.18 of the total. The late epoch 
declines are also different, the rate scaling as and G®-* for 
the DDC and DDW models respectively. 

Fig. 1^ shows the predicted rates as a function of redshift 
for each of the three models, having assumed the Madau and 
Dickinson cosmic SFH. The best ht of the models with obser¬ 
vations was derived by least square minimisation weighting the 
measurement by their cr and gives Ki^ - 7.5 x 10“"^ M"' for SD 
and DDW and 8.5 x lO^^ Mg' for DDC. 

The smallest residuals in the whole redshift range is ob¬ 
tained for the SD model (rms= 0.0028). The DDC model (rms= 
0.0088) gives an excellent fit up to redshift z ~ 1.2 but predicts a 
relatively mild decline at higher redshift (still not in conflict with 
the observations). Instead, the DDW model (rms= 0.012) shows 
an overall shallow evolution with respect to the observations. 
Our conclusion is that since the dispersion of the rate measure¬ 
ments is comparable to the scatter of theoretical tracks we cannot 
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discriminate between SD and DD models, tough for the DD sce¬ 
narios model with close binary separation seems favoured. 

We remark that the impact of adopting different cosmic 
SFHs, MD14 or HB06, on the predicted rate is negligible, while 
it is more relevant for the estimate of constant Ki^, that is the 
number of stars which end up as a SNIa per unit mass of the 
parent stellar population. In particular, fitting the observations 
with HB06 SFH requires = 5.9, 5.7, 6.9 xlO^^ for SD, 
DDW and DDC respectively, values which are ~ 20% smaller 
than those obtained when using the MD14 SFH. The number of 
potential progenitors per unit stellar mass depends on the IMF 
and, assuming from to 8 Mq range of SN la progenitors, it is 
0.021 for a plain Salpeter IMF and 0.028 for a SalA IMF. Stars 
in the selected mass range that should end up as SN la t oaccount 
for the observed rates are 4% and 2% using the MD14 and HB06 
SFH, respectively. These fractions are close to the lower edge of 
the range reported in Maoz and Mannucci (2012). While they 
remain large with respect to most theoretical predictions from 
binary population synthesis codes, they still represent a minor 
fraction of all potential progenitors. 

In any case, we conclude that there is no need to invoke ad- 
hoc delay time distributions, unrelated to the standard expecta- 
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redshift 


Fig. 20. Average values of the SN la rates per unit volume as a func¬ 
tion of redshift. The still unique measurement at z > 2 of [Rodney et al.| 
( |2014^ is also plotted. For the derivation of the SNIa rate model evolu¬ 
tion we adopted the SFH from|Madau & Dickinson|(|2014|l. 


tions from stellar evolution theory. This confirms earlier conclu¬ 
sions (e^^gjF&steFeFar|20^ Blanc & GreggiopOOS Greggio 


let al.|2008t. 


In a forthcoming paper (PII) we will use a detailed charac¬ 
terization of the properties of the galaxy sample in order to in¬ 
vestigate the dependence of SN rates on galaxy parameters, and 
to obtain additional constraints on the SN progenitor scenarios. 


9. Conclusions 


We presented the preliminary results of a new SN search, SU¬ 
DARE, that was designed to measure SN rates in the redshift 
range 0 < z < 1. This paper describes the survey strategy, the 
selection and confirmation of candidates, the construction of the 
galaxy catalog and the rate estimates based on the first two years 
of the survey. 

As search fields we selected two of the best studied extra- 
galactic fields, namely CDFS and COSMOS for which a wealth 
of multi-band coverage is available. Our own data, the syn¬ 
ergy with the VOICE project, the complementary data from the 
VIDEO survey, along with public data from the literature al¬ 
lowed us to obtain a multi-band photometric catalog for the 
galaxies with magnitudes K < 23.5 that was exploited to esti¬ 


mate the photometric redshift using the EAZY code (Brammer 
let al.|2008» . 

We discovered 117 SNe, of which 27 were assigned a 
weight= 0.5 due to a poor template match or to a low number 
of detections. Most of the SNe are classified of type la (57%). 
Eor the core collapse, 44% are type II, 22% type Iln and 34% 
type Ib/c. 

With this SN sample and an accurate measurement of the de¬ 
tection efficiency of our search, we computed the rate of SNe 
per unit volume. Eor the CC SNe our measurements are in 
excellent agreement with previous results and fully consistent 
with the predictions from the cosmic SFH of Madau & Dickin- 
|son| ( |2014| , assuming a standard mass range for the progenitors 
(8 < M < 40 Mq). Therefore, previous claims of a significant 
disagreement between SFH and CC SN rates are not confirmed. 
This conclusion relies on the revision of the cosmic SFH be¬ 


cause our rate estimates are consistent withother measurements 
from literature. 

For the SN la, our measurements are consistent with liter¬ 
ature values within the errors. We conclude that the dispersion 
of SN la rate estimates and the marginal differences for the evo¬ 
lution with cosmic time of the volumetric SN rate does not al¬ 
low us to discriminate between SD and DD progenitor scenarios. 
However, with respect to the three tested models(SD, DDC and 
DDW from Greggio|2010) l, the SD gives a better fit on the whole 
redshift range whereas the DDC appears to perfectly match the 
steep rise of the rate up to redshift 1.2. The DDW model that cor¬ 
responds to a wide binary separation and a relatively flat delay 
time distribution appears disfavoured. 

As a first attempt at searching for evolution of SN diversity, 
we found no evidence of evolution of the SN Ib/c fraction. The 
fraction of type Iln SNe detected in the 0.15 < z < .35 redshift 
bin is consistent with the measurements for the local Universe. 
The rate in the higher redshift bin is formally significantly lower. 
Whether this is evidence for some evolution or a bias in our sur¬ 
vey needs to be verified with more data. 
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